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Abstract The aim of this paper is to use large deviation theory in order to compute the entropy of 
macrostates for the microcanonical measure of the shallow water system. The main prediction of this 
full statistical mechanics computation is the energy partition between a large scale vortical flow and small 
scale fluctuations related to inertia-gravity waves. We introduce for that purpose a semi-Lagrangian dis¬ 
crete model of the continuous shallow water system, and compute the corresponding statistical equilibria. 
We argue that microcanonical equilibrium states of the discrete model in the continuous limit are equilib¬ 
rium states of the actual shallow water system. 

We show that the presence of small scale fluctuations selects a subclass of equilibria among the states 
that were previously computed by phenomenological approaches that were neglecting such fluctuations. In 
the limit of weak height fluctuations, the equilibrium state can be interpreted as two subsystems in thermal 
contact: one subsystem corresponds to the large scale vortical flow, the other subsystem corresponds to 
small scale height and velocity fluctuations. It is shown that either a non-zero circulation or rotation and 
bottom topography are required to sustain a non-zero large scale flow at equilibrium. 

Explicit computation of the equilibria and their energy partition is presented in the quasi-geostrophic 
limit for the energy-enstrophy ensemble. The possible role of small scale dissipation and shocks is dis¬ 
cussed. A geophysical application to the Zapiola anticyclone is presented. 
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I Introduction 

Geophysical turbulent flows have the propensity to self-organize into large scale coherent structures such 
as cyclones, anticyclones and jets. These coherent structures are long lived, but can also loose energy, 
for instance through the radiation of waves that eventually break into small scale structures. The aim of 
this paper is to understand the energy partition into large scale structures and small scale fluctuations in 
the framework of freely evolving shallow water dynamics, using statistical mechanics arguments. Indeed, 
geophysical turbulent flows involve a huge number of degrees of freedom coupled through non-linear 
interactions, which strongly motivates a statistical mechanics approach. This approach allows to reduce 
the study of self-organization and energy partition down to a few parameters, such as the total energy of 
the flow and its total enstrophy. 

In the case of the three dimensional Euler equations, equilibrium statistical mechanics predicts that all 
the energy is lost into small scales, consistently with the classical picture of a small scale energy transfer. 
By contrast, two-dimensional flows are characterized by a large scale energy transfer, and equilibrium 
tools are appropriate to describe the large scale structure resulting from self-organization at the domain 
scale, in the absence of forcing and dissipation. The idea goes back to Onsager ll28l . and has been mostly 
developed during the nineties after the work of Miller-Robert-Sommeria 112311331 . see also Refs. 111711 101 

II and references therein. Importantly, the theory predicts that the contribution of small scale fluctuations 
to the total energy are negligible in the two-dimensional case. Equilibrium statistical mechanics of two- 
dimensional and quasi-geostrophic flows is now fairly well understood. It has been applied to various 
problems in geophysical context such as the description of Jovian vortices 1144 l l3l. oceanic rings and jets 
B501I55I . equilibria on a sphere lfl5l , and to describe the vertical energy partition in continuously stratified 
quasi-geostrophic flows IU9ll52ll48l . 

Due to the combined effect of stable stratification, thin aspect ratio and rotation, geophysical flows 
are very different from classical three-dimensional turbulence. However, such flows are not purely two- 
dimensional. Here we consider the shallow water equations, which is an intermediate model between 
three-dimensional and two-dimensional turbulence. This model describes the dynamics of a thin layer 
a fluid with homogeneous density. On the one hand, shallow water equations admit conservation laws 
similar to two-dimensional Euler equations, that lead to self-organization of the energy at large scale in 
the Euler case. On the other hand, shallow water dynamics support the presence of inertia-gravity waves 
that are absent from purely two-dimensional turbulence. A small scale energy transfer may exist due to 
the existence of these inertia-gravity waves in the shallow water system. The quasi-geostrophic model is 
recovered as a limit case of the shallow water model, when the Rossby parameter (comparing inertial terms 
to Coriolis forces) is small. A small Rossby number corresponds to a strong rotation limit. It is then natural 
to ask whether previously computed statistical equilibrium states of the quasi-geostrophic models remain 
equilibrium states of the shallow water model. More generally, given a certain amount of energy in an 
unforced, undissipated geophysical flow, will the flow self-organize into a large scale coherent structure, 
just as in two dimensional turbulence? Or will the energy be transferred towards the small scales, just 
as in three dimensional homogeneous turbulence? The aim of this paper is to answer these questions by 
computing statistical equilibrium states of the inviscid shallow water model. 

The first step before computing equilibrium states is to identify a suitable phase space to describe 
microscopic configurations of the system. The phase space variables must satisfy a Liouville theorem, 
which ensures that the flow in phase space is divergence-less. Consequently, a uniform measure on a 
constant energy-Casimirs shell of phase space is invariant (microcanonical measure). The second step is 
to describe the system at a macroscopic level. The macrostates will be the sets of microstates sharing 
the same macroscopic behavior. The third step is to find the most probable macrostate, and to show that 
almost all the microstates correspond to this macrostate for given values of the constraints. While these 
three steps may be straightforward for differential equations with a finite number of degrees of freedom, for 
continuous systems described by partial differential equations, these three steps require the introduction of 
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discrete approximations of the continuous field and of the invariant measure, and to study the continuous 
field limit of these discrete approximations. This point will be further discussed in the following. 

This program has been achieved in the 90s for the two-dimensional Euler equations. Indeed, a Liouville 
theorem is satisfied by the vorticity field, which describes therefore a microscopic configuration of the 
system. A macrostate can be defined as a probability field describing the distribution of vorticity levels at 
a given point, either through a coarse graining procedure 12211231 or directly by the introduction of Young 
measures 00113111331 . 

For the shallow water dynamics, they are further issues that need to be overcome. The existence of a 
Liouville theorem for the shallow water flow was found by Warn l!53l by describing the flow configurations 
on a basis given by the eigenmodes of the linearized dynamics. However, the constraints of the problem 
given by dynamical invariants are not easily expressed in terms of the variables satisfying this Liouville 
theorem (except in the weak flow limit discussed by Warn [[53]]). This difficulty has been overcome by 
Weichman and Petrich 1561 who considered first a Lagrangian representation of the flow, and then used a 
formal change of variable to describe the flow configurations with Eulerian variables convenient to express 
the constraints of the problem. Using a different method that does not require a Lagrangian representation 
of the fluid, we will show the existence of a formal Liouville theorem. 

A second difficulty concerns the choice of a relevant discrete approximation that allows to keep as 
much as possible geometric conservations laws of the continuous dynamics. Those geometric conserva¬ 
tion laws include the Liouville property, the Lagrangian conservation laws (i.e. the conservation of the 
volume carried by each fluid particle), and global dynamical invariants. Unfortunately, we are not aware 
of a discrete model that does not break at least one of those geometric conservation laws. However, we 
will argue that there is no logical need for the discrete model to satisfy all the conservation laws of the 
continuous dynamics in order to guess the correct microcanonical measure of the continuous system by 
considering the limit of a large number of degrees of freedom. 

A third difficulty is that local small scale fluctuations of the fields may have a substantial contribution to 
the total energy. This contrasts with 2d Euler dynamics, for which small scale fluctuations of the vorticity 
field have a vanishingly small contribution to the total energy in the continuous limit. In the shallow water 
case, it is not a priori obvious that the contribution of these small scale fluctuations to the total energy can 
be expressed in terms of the macroscopic probability density field. 

In order to overcome the second and the third difficulty, we introduce in this paper a semi-Lagrangian 
discretization : one the one hand, we consider fluid particles of equal volume, which allows to keep track of 
the Lagrangian conservation laws of the dynamics. On the other hand, the particle positions are restricted 
to a uniform horizontal grid, and each grid node may contain many particles, which allows for an Eulerian 
representation of the macrostates and to keep track of global conservation laws, while taking into account 
the presence of small scale fluctuations contributing to the total energy. We then derive the statistical 
mechanics theory for this discrete representation of the shallow water model, using large deviation theory. 

We will argue that in the limit of a large number of degrees of freedom, the equilibrium states of the 
semi-Lagrangian discrete model are the equilibrium states of the actual shallow water system, and we 
will also show that we recover with this model results already obtained in several limiting cases. While 
the statistical mechanics treatment of our discrete model is rigorous, there is some arbitrariness in our 
definition of the discrete model, which is not fully satisfactory. 

We stress that our approach to guess the invariant measure of the continuous shallow-water system is 
heuristic: there is to our knowledge no simple way to define mathematically the microcanonical measure 
or the Gibbs measure of an Hamiltonian infinite dimensional system. As far as we know the only example 
of a rigorous work for defining invariant measures in the class of deterministic partial differential equations 
of interest, is a work by Q on the periodic nonlinear Schrodinger equation. For the 2d Euler equations, 
the microcanonical measure seems clear from a physical point of view because different discretizations, 
for instance either Eulerian ones B22U321I or Lagrangian ones (see (81, lead to a consistent picture in the 
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thermodynamics limit, but even in that case no clear mathematical construction of the invariant measure 
exists. 

There exist only a very limited number of results on statistical equilibrium states of the shallow water 
system. Warn (53l studied the equilibrium states in a weak flow limit in the energy-enstrophy ensemble. 
He showed that in the absence of lateral boundaries and in the absence of bottom topography, all the 
initial energy of a shallow water flow is transferred towards small scales. Here we relax the hypothesis of 
a weak flow, generalizing the conclusions of Warn 15311 for any flow, and we discuss the effect of bottom 
topography. We show that when there is a non zero bottom topography and when the flow is rotating, there 
is a large scale flow associated with the equilibrium states. 

Refs. 1201 171 did compute statistical equilibrium states of "balanced” shallow water flows, by assuming 
that all the energy remain in the large scale flow. The equilibria described by Merryfield et al l20l were ob¬ 
tained in the framework of the energy-enstrophy theory, neglecting any other potential vorticity moments 
than the potential enstrophy. Similar states were described as minima of the potential enstrophy for the 
macroscopic potential vorticity field by Sanson l37l . 

Chavanis and Sommeria 0 proposed a generalization of the 2D-Euler variational problem given by the 
Miller-Robert-Sommeria theory to the shallow water case. Their main result is a relationship between the 
large scale streamfunction and height field. This result was very interesting and inspiring to us. However, 
Chavanis and Sommeria did not derive their variational problem from statistical mechanics arguments but 
proceeded through analogy. They were moreover neglecting height fluctuations. This does not allow for 
energy partition between vortical flow and fluctuations, and moreover leads to inconsistencies for some 
range of parameters, as the Chavanis-Sommeria constrained entropy can be shown to have no maxima 
for negative temperatures (the negative temperature critical points of this variational problem are saddles 
rather than maxima). As a consequence a range of possible energy values can not be achieved in this 
phenomenological framework. Nevertheless, with an approach based on statistical mechanics, we will 
confirm, in this study, the form of the variational problem proposed by Chavanis and Sommeria (7) for 
describing part of the field, for a restricted range of parameters. We stress however, that generically part of 
the energy will be carried by the fluctuations, and thus that the Chavanis-Sommeria variational problem 
should be considered with an energy value which is not the total one. Determining which part of the 
energy should be taken into account requires a full statistical mechanics treatment, taking into account 
height fluctuations. We also note that for some other ranges of parameters, all the energy will be carried 
by the fluctuations, and thus the Chavanis-Sommeria variational problem then does not make sense. 

In the preparation of this work, we have also been inspired by the work of Weichman and Petrich 1561 . 
In this work, the authors computed a class of statistical equilibrium states of the shallow water system, 
starting from a grand canonical distribution. The main result of their work is the same equation describ¬ 
ing the relation between the large scale stream function and height field, as the one previously obtained by 
Chavanis and Sommeria 0 . The Weichman-Petrich approach seemed more promising than the Chavanis- 
Sommeria one as it was based on a statistical mechanics treatment. However it also fails to predict which 
part of the energy goes into fluctuations, and to recognize the range of parameters for which the mean field 
equation for the largest scale is relevant. The reason is that, while natural in condensed matter physics, 
the hypothesis of a Gibbs or a grand canonical distribution, is questionable for an inertial flow equation 
which is in contact neither with an energy bath, nor with a potential vorticity bath. Related to this issue, 
Weichman and Petrich had to assume an ad hoc scaling for the thermodynamical parameters, in order to 
obtain statistical ensembles where entropy actually balances conservation laws. It has been recognized for 
a long time that this kind of theoretical difficulties are related to the Rayleigh-Jeans paradox. The proper 
way to address these issues, in some fluid models like for instance the two-dimensional Euler equation, is 
to start from the microcanonical measure rather than from a canonical or grand canonical one (41. We will 
overcome all these problems in this work by solving the microcanonical problem, which may seem more 
difficult, but which can be handled using large deviation theory. It also leads to a more precise description 
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of macrostate probabilities. 

In this paper, starting from the microcanonicai measure, we propose a complete computation of the 
macrostate entropy for the shallow water equations. In order to achieve this goal, we consider a large- 
deviation approach that allows to perform the statistical mechanics computation in an explicit and clear 
way, and that gives a precise description of most-probable macrostates. Moreover, we compute explicitly 
both the large scale flow and the small scales fluctuations of the equilibrium states. Our results are thus a 
complete statistical mechanics treatment of the shallow water equations. 

We also connect our general results to some of the partial treatment discussed in the previous litera¬ 
ture. For instance, we show that only one subclass of the states described by 112011561171 are actual statistical 
equilibria of the shallow water model. More precisely, we will see that the large scale flow and the small 
scales fluctuations can be interpreted as two subsystems in thermal contact, and that the temperature of 
small scales fluctuations is necessarily positive. The large scale flow is therefore also characterized by a 
positive temperature. 

The paper is organized as follows. The shallow water model and its properties are introduced in section 
El Equilibrium statistical mechanics of a discrete flow model arguably relevant to describe the continuous 
shallow water system is derived in section [3] Computation of the equilibrium states and a description of 
their main properties are presented in section [4] Energy partition between a large scale flow and small 
scale fluctuations is computed analytically in the quasi-geostrophic limit in section[3 which also includes 
a geophysical application to the Zapiola anticyclone. The main results are summarized and discussed in 
the conclusion. 


2 Shallow water model 


2.1 Dynamics 


The shallow water equations describe the dynamics of a fluid layer with uniform mass density, in the 
limit where the layer depth is very small compared to the horizontal length scales of the flow B461I29I . In 
this limit the vertical momentum equation yields hydrostatic equilibrium, and the horizontal velocity field 
(u, v) = u(x,t) is depth independent, where x = (x,y) is any point of a two-dimensional simply connected 
domain £0. We consider the Coriolis force in the f-plane approximation, i.e. with a constant Coriolis 
parameter / and a rotation axis along the vertical direction. We denote t] (x, t) the vertical displacement of 
the upper interface and hf, (x) the bottom topography (see Fig.|T). The origin of the z axis is chosen such 
that 

f dx h b (x) = 0, (1) 

J® 

and the vertical displacement is defined such that 


/ 

7 ® 


dx r] (x,f) = 0 


as well. We introduce the total depth 

h—H-hb + r] (x,t), 

where 



( 2 ) 

(3) 

(4) 
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Fig. 1: Scheme of a vertical slice of fluid, T) is the upper interface displacement, /;/, is the bottom topogra¬ 
phy, h is the height field and H is the mean height of the fluid ((h)) 


is the mean depth of the fluid, with \S>\ the area of the flow domain. The horizontal and vertical length 
units can always be chosen such that the domain area and the mean height H are equal to one (\S>\ — 1. 
H = 1), and this choice will be made in the remainder of this paper. 

The dynamics is given by the horizontal momentum equations 


d t u + ud x u + vdyU — fv — —gd x (h + hb) 
d t v + ud x v + vdyV + fu — —gd y (h + h h ) 


and the mass continuity equation 


(5) 


d,h + V • (hi l) = 0, 

where V = (d x ,d y ), with impermeability boundary conditions 

u n = 0 on dSt, 

where n is the outward border-normal unit vector. 

By introducing the Bernoulli function 

B= ^u 2 +g(h + h b - 1), 


( 6 ) 

(7) 

( 8 ) 




the relative vorticity field 


(9) 


co = d x v — d y u, 


the divergence field 


£ = V • u = d x u + d y v, 


and the potential vorticity field 


q = 


<» + / 

h 


the momentum equations (0 can be recast into a conservative form 


( 10 ) 

(ID 


c),u qhv — —d x B 
r),v + qhu — —d y B 


( 12 ) 


One can get from the momentum equations (1121) a dynamical equations for the potential vorticity field (II 11 1 
and for the divergence field ( fTOl ). respectively: 


d t q + u'Vq — 0, 


03) 


e>,£ - V x - (>/m) = -AB. (14) 

where A = V • V is the Laplacian operator and V = (— d y ,d x ). In order to show that the shallow wa¬ 
ter system described by (/i,u) is also fully described by the triplet (h.q. £), we introduce the velocity 
streamfunction \jf(x,y,t) and the velocity potential cj) (x, y. I ) through the Helmholtz decomposition of the 

velocityQ 


u = V J -y/' +V0. 

The boundary condition for the velocity given in Eq. (0) yields 


n 'V J -\j/ = 0 
n • V0 = 0 


on d@, 


(15) 


(16) 


where n is the boundary normal vector. The relative vorticity © and the divergence (1101) can be expressed 
in terms of the velocity streamfunction and velocity potential as 


C0 = A\f/ 7 ^ —A(j) . (17) 

To the fields co and £ correspond a unique field if/ and a field <j> defined up to a constant with the boundary 
condition (IT6l) . Thus, to the fields CO = qli — f and £ corresponds a unique velocity vector written formally 
as 

u = V- l 4 -1 ($/i-/)+V 4 _i £. (18) 

It will also be useful to consider the field u defined as 

B = A l/2 <\>, (19) 

where A is the linear operator whose eigenmodes are Laplacian eigenmodes in the domain c i/) with the 
boundary condition n • V = 0 on dS>, and whose eigenvalues are the negative square root of the modulus 

1 In the particular case of a bi-periodic domain, i.e. with periodic boundary conditions for the velocity u, one would need 
to describe in addition the homogeneous part of the velocity field, which is both divergence-less and irrotational. 
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of Laplacian eigenvalues. The field fi can be interpreted as a measure of the amplitude of the potential 
contribution to the velocity field, given that f^ju 2 — f@(V<t>) ■ This field will also be referred to as the 
divergence field in the following. To the field /u correspond a unique field (j) (up to a constant). Thus, 
the shallow water system is fully described by the triplet (,), and the velocity vector can be written 
formally in terms of these fields as 


u = V ± A~ 1 [hq -/] + WV- :1 / 2 [jU]. 


( 20 ) 


2.2 Stationary states 


We investigate in this subsection the conditions for stationarity of the flow in the shallow water model fol¬ 
lowing (7J. Using a Helmholtz decomposition, one can define the mass transport streamfunction 0(x,y,t) 
and the mass transport potential <P(x,y,t) , not to be confused with \ff and (!) defined in Eq. (ITHT i and (fl6l> : 


hu — V x 0 +V0. 

The boundary condition for the velocity given in Eq. (|7|) yields 


n • V x 0 = 0 
n V0 = 0 


on dS). 


( 21 ) 


( 22 ) 


It will be useful in the remainder of the paper to express the vorticity field defined in Eq. ([9} and the 
divergent field defined in Eq. dTTTb in terms of these transport streamfunction 0 and transport potential <P : 


<P 


C = v (^)-./(~0 


where 


J (/, g) = djdyg - d y fd x g 


(23) 

(24) 

(25) 


is the Jacobian operator. The dynamics in Eqs. ©, (fl3l) and ( fl4l l can also be written in terms of W and 0: 

d,h — —A<P, (26) 


d t q+^J(W,q) + ^V0-Vq = 0, (27) 

d r £+J(cp,q)-qA'F-Vq-W = -AB. (28) 

We see that d,h — 0 implies 0 = Cst, and d t q — 0 implies J(W,q) = 0, which means that isolines of 
potential vorticity are the mass transport streamlines. This is the case if for instance q — F (¥). Recip¬ 
rocally, if 0 = Cst and if q and 0 have the same isolines, then d,h — 0 and d,q = 0, which also implies 
d,W — 0. Thus, through the decomposition ( fTil . the velocity field is also stationary, with u = (I fh) V 1 0. 
We conclude that a necessary and sufficient condition for a shallow water flow (h, u) to be stationary is 


0 = Cst, q and 0 have the same isolines. 


(29) 
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There is an additional relation verified by the stationary flow. This relation may be obtained by consid¬ 
ering stationarity of the kinetic energy. The dynamics of the kinetic energy is obtained from Eq. (112b : 





--J{W,B)~-V<P'VB. 
h h 


(30) 


Stationarity of the kinetic energy field and Eq. ( |291 ) gives 7( t / / . B) = 0: the Bernoulli function defined in 
Eq. © and the mass transport streamfunction W have the same isolines. In addition, in any subdomain 
where q = F(W), the stationarity of the velocity field and the momentum equations (112b give the relation 


q = 


AB 

d*F' 


(31) 


2.3 Conserved quantities 


Provided that the velocity and height fields remain differentiable, the shallow water dynamics in Eqs. ©. 
( [T3T > and (THT i conserves the total energy 


£ [u,/t] = ■ 


dx 


lm 2 +g(h + h b - 1)" 


(32) 


which includes a kinetic energy contribution and a potential energy contribution. It is known that the shal¬ 
low water dynamics sometimes leads to shocks that prevent energy conservation. We postpone a discussion 
of this important point to the last section of this paper. 


As a consequence of particle relabeling symmetry there is an infinite number of other conserved quan¬ 
tities called Casimir functionals, see e.g. ED- These functionals are written 


% = j dx hs(q ) , (33) 

where 5 is any distribution. This conservation can be easily checked from Eqs. (113b and ©. The conserva¬ 
tion of all the Casimirs implies the conservation of all the potential vorticity moments 


VJt e N, 2k 


dx hq k . 


(34) 


These Casimir functionals include the total mass conservation (k — 0), the conservation of the circulatiorQ 
(k — 1) and the conservation of the enstrophy (k — 2). 

The conservation of all the Casimirs is equivalent to the conservation of mass plus the conservation of 
the potential vorticity distribution @ [q\ defined through 


Vcr, S>[q]{o)Ao ^ j Axh^ {a < q < a+Ao} (35) 

where J^{ a <q<a+da} i s the characteristic function, i.e. it returns one if o <q(x) < O + d<7 and zero 
otherwise. It means that the global volume of each potential vorticity level a q is conserved through the 
dynamics. 

2 The actual circulation is usually defined as F = f 3D dl • u, where dl is a vector tangent to the domain boundary. Stokes 
theorem yields to Z\ = F + /. 
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We will restrict ourselves to initial states where the potential vorticity field has a distribution charac¬ 
terized by its moments only, then the knowledge of its global distribution given in Eq. (1351) and of the total 
mass is equivalent to the knowledge of the moments given in Eq. (1341 ). 

Depending on the domain geometry, there could be additional invariants. For the sake of simplicity, we 
do not discuss the role of these additional invariant in this paper, but it would not be difficult to generalize 
our results by taking them into account. 

In the quasi-geostrophic or two-dimensional Euler models, dynamical invariants have important conse¬ 
quences such as the large scale energy transfer, the self-organization at the domain scale, and the existence 
of an infinite number of stable states for the dynamics, see for instance Ref. f4j| and references therein. 
We show in the present paper that these dynamical invariants play a similar role in the shallow water case, 
allowing for a large-scale circulation associated with the potential vorticity field, even if this process may 
be associated with a concomitant loss of energy toward small scales. 


3 Equilibrium statistical mechanics of a discrete shallow water model 

The aim of this section is to compute statistical equilibrium states of a discrete model of the shallow 
water system and to consider the thermodynamic limit for this model. In the first subsection, a formal 
Liouville theorem is given for the triplet of variables ( h,hu,hv) and then the triplet (h.q. ii) after a change 
of variables. The derivation given with more details in Appendix |A]is made in the Eulerian representation. 
This allows to write formally the microcanonical measure of the shallow water model for these triplets of 
fields. 

In the second subsection, a finite dimensional semi-Lagrangian discretization of the model is proposed 
to give a physical meaning to the formal measure. Other discretizations of the shallow water system may 
have been considered but any other choices we have tried were not taking into account the fluid particle 
mass conservation and led to inconsistent results (for instance the equilibrium states were not stationary, 
not stable by coarse-graining, and unbalanced, see Appendix [Kj. 

Then the macroscopic states are defined through a coarse-graining procedure in a third subsection. 
The variational problem leading to the equilibrium states, i.e. the most probable macroscopic state, is 
introduced in the fourth subsection. This variational problem generalizes the Miller-Robert-Sommeria 
equilibrium theory to the shallow water system. 


3.1 Liouville theorem 

The first step before computing equilibrium states is to define what is a microscopic configuration of the 
system, which requires to identify the relevant phase space. The simplest set of variables to consider are 
those that satisfy a Liouville theorem. A Liouville theorem means that the flow in phase space is non- 
divergent, which implies the time invariance of the microcanonical ensemble. The existence of a Liouville 
theorem for the shallow water system was initially shown by Warn lf53l who considered a decomposition 
of the flow fields on a basis given by the eigenmodes of the linearized dynamics. A Liouville theorem 
was shown by Weichman and Petrich ll56l for a Lagrangian representation of the dynamics. Following a 
general method proposed by one of us and described by Thalabard l40l , it is shown in Appendix[A]that the 
triplet (h,hu,hv) does satisfy a formal Liouville theorem. At a formal level, the microcanonical measure 
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can then be writter@ 


d HhfiuM (E,{Zk} k > 0 ) = ^ [h] ® M 9 [hv] 8 {S - E) fj 5 (J* - Z k ), (36) 

Ll [h.,{L k \ k >Q) k=Q 

with the phase space volume 

Q (E, (Z,},> 0 ) = f 9 [h] 9 [hu] $> [hv] 5 (<?-£) f] 5 (S’* - Z*). (37) 

■I k=0 

Here S is the energy of a microscopic configuration defined in Eq. (1321) . and the {3f k } k> o are the potential 
vorticity moments of a microscopic configuration defined in Eq. (134b . Those constraints are the dynamical 
invariants of the shallow water model. The notation f 9 [h] 9 [hu] 9 [hv] means that the integral is formally 
performed over each possible triplet of fields (h, hu, hv). The microcanonical measure allows to compute 
the expectation of an observable A [hjiu. hv] in the microcanonical ensemble as 


^duhMM = / d Vh,huM */[h,hu,hv ]. 


(38) 


Assuming ergodicity, the ensemble average (sZ) d ^ can finally be interpreted as the time average of 
the observable si. 


The triplet ( h,hu,hv ) is not a convenient one to work with, since the Casimir functionals {'Z'k} k >() 
defined in Eq. ( 134b are not easily expressed in terms of these fields. Indeed, the expression of the Casimir 
functionals {2f k } k >o involve not only the triplet (h.hu.hv), but also the triplet of the horizontal derivatives 
of these fields. We showed in subsection 12. II that the triplet of fields (h,q, jU) fully describes the shallow 
water dynamics in a closed domain, and the functionals {3f k } k >o are much more easily expressed in terms 
of the fields h,q. It is therefore more convenient to use these fields as independent variables. The price 
to pay is that the simple form of the energy defined in Eq. (f32l> in terms of the triplet (li.hu. vh) becomes 
more complicated when expressed in terms of the triplet (li,q,p.). However, we will propose a simplified 
version of this energy functional, and argue in Appendix [B] that this is the relevant form of the energy to 
consider to compute the equilibrium state. 


Unfortunately no direct proof of a Liouville theorem can be obtained for the triplet (h,q,/J,). However, 
it is still possible to start from the microcanonical measure built with (h,hu,hv) in Eq. ( 1361 ). and change 
variables at a formal level.lt is shown in Appendix I A. 2 1 that the Jacobian of the transformation isQ 


J 


(h,liu,hv) 

. (q,h,n) . 


— Civ 1 . 


(39) 


where C is a constant. We conclude that the microcanonical measure can therefore be formally written 


d Hh.q.ji {E, {Z k } k > 0 ) = 1 J i 3 4 ® [h] S [q] 9 [M] 5 {g - E) ff 8 (3T k - Z k ), (40) 

Ll (£,; A=0 

3 The letter “p” appearing in the measure denoted dp is not related to the divergent field denoted p. 

4 The fi 3 term that appears after the change of variables in the functional integral must be understood as a “functional 
product” nxesTi 3 (x), see also the finite-dimensional representation of this measure given in the next subsection. 













13 


with 

<2(£,{Z*}*> 0 )= f h^[h\®[q}®[v}5{£-E)Y\5{2? k -Z k ) . (41) 

J k =0 

In that case, the expectation of an observable A [q 7 h,p] in the microcanonical ensemble is given by 

^dn qM = j d n qAll */[q,h,H] ■ (42) 

A finite-dimensional projection of the fields will be given in the next subsection to give a meaning to the 
formal notations of this subsection. 


3.2 A discrete model 

In this subsection, we devise a discrete model of the shallow water system based on a semi-Lagrangian 
representation. This allows to give a finite dimensional representation of the formal measure given Eq. 
®. 

In order to keep track of the conservation properties of the continuous dynamics, we propose a con¬ 
venient discretization of the shallow water system in terms of fluid particles. Since the fluid is considered 
incompressible, it is discretized into equal volume particles. The fluid particles do not have necessarily the 
same height, and the horizontal velocity field may be divergent. It is not possible to build a uniform grid 
with a single fluid particle per grid point since two particles with equal volume and different height can 
not occupy the same area. In order to bypass this difficulty, we define a uniform grid where several fluid 
particles can occupy a given site. To keep track that the actual fluid is continuous and incompressible, we 
add the condition that the area occupied by the particles inside a grid site fit the available area of the grid 
site. 

In a first step the ensemble of microscopic configurations (the ensemble of the microstates of the 
discrete model) is defined. In a second step, the constraints of the discrete shallow water system are intro¬ 
duced, which allows to define the microcanonical measure of the discrete model in a third step. 


3.2.1 Definition of the ensemble of microscopic configurations 


For the sake of simplicity, the domain S> where the flow takes place is considered rectangular but gener¬ 
alizing the results to any shape would be straightforward. We recall that the horizontal and vertical length 
units have been chosen such that the domain area and the mean height are equal to one (\2>\— L x L y — 1, 
H = 1). The domain S> is discretized into a uniform grid with N = N x x N y sites. The area of a grid site 
is \S>\ /N = 1 /N. Each site can contain many fluid particles. It is assumed that the fluid contains N x M 
fluid particles of equal mass and volume (the fluid is incompressible), where M is the average number of 
particles per site. The volume of a fluid particle is therefore 


8v = \3E = J_. 

NM NM 


(43) 


The grid sites are indexed by (i, j) with 1 < i,j < N x .N y and the fluid particles by n with 1 < n < NM. 


Each fluid particle is characterized by its position (/„,/„) on the grid, by its potential vorticity q n 6 
[q min, q max], its divergence ,U„ <E and its height h n e [h min .h max \. The cutoffs on the potential 

vorticity can be physically related to actual minimum and maximum in the global distribution of potential 
vorticity level defined in Eq. ( 1351 ). since this distribution is conserved by the dynamics. Such a justification 
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does not exist for the other cutoffs. We will consider first the limit of an infinite number of fluid particles 
per grid site (M —> +°°), then the limit of an infinite number of grid site (N —> <*>) and finally the limit of 
infinite height and divergence cut-off jl max —> +°°, h max —> +°°, -4 0. We will see that the result does 

not depend on those cut-off. 


Let us introduce Mu the number of particles per grid site (i,j), defined as 

NM 

Mi, - I SlnAnJ- (44) 

n= 1 


The set of the particles that belong to the site (i, j) is denoted 


J(ij — {1 < n <NM\ = ( i,j )}, 


(45) 


whose cardinal is Mjj. Mass conservation states that the total number of particles filling the grid is a 
constant, which gives the constraint 

Y / M iJ =NM, (46) 

ij 


where means that we sum over the all the sites of the grid. A fluid particle labeled by n and carrying 
the height h n occupies an area SV jh n . The constraint that the area of each grid site is covered by fluid 
particles leads to the constraint 


Vi, 7 


- y -i = i, 

M ^ h 




.. hn 


(47) 


where .Mij is the set defined in Eq. (|45T >. 

The ensemble of microstates of the discrete model is given by the set of all reachable values of grid po¬ 
sitions, potential vorticities, divergences and heights of each fluid particle in accordance with the constraint 
of particle filling the area of each grid site: 


Xmicro — A Xmicro — {(//;, 7/, ) ■ (ijn ■ f-i-n ■ h n j | < n < ^ y 


Vi, 7 1 < i,j<N x ,N y 


1 

M 




\. hn 


= 1 


(48) 


3.2.2 Coarse-graining 

Here we consider a microstate Xmicro — {(JnMn) ■ <-ln■ Hn ■ I'n} i m n< MN ani -( an arbitrary function 

8 '■ n—¥g n = g{h n ,q n ,n„). (49) 

We introduce two different coarse-graining procedures: an areal coarse-graining, and a volumetric coarse- 
graining. 

The areal coarse-graining of the function g is defined at each grid point (i. j) as 


_ 1 y 1 

~ M h Sn 




hn 


(50) 


where .cXX-u is a set defined in Eq. (145 b . The terms l/h n appearing in Eq. (150b means that we consider 
local average of g n weighted by the area occupied by each fluid particle. Note that we will only consider 
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function g such that ~g tJ converges to a finite value in the limit of large Mjj ~ M. This means that the terms 
g n should not be allowed to scale with M. This is the reason why we will consider first the large M limit, 
and then the limit of large cut-off jl max and h max for the fields Li and h. 

The area filling constraint in Eq. (147b can then be written in terms of this areal coarse graining: 


1 ij — 1 ■ (51) 

for any i,j. We also notice that the areal coarse-grained height field is simply the ratio of the number of 
particles in the site (i,j) over the averaged number of particles per site: 


hij 


Mil 

M ' 


(52) 


The volumetric coarse-graining of the function g is defined at each grid point (i, j) as 


<*>:y = W E 


gn- 


l J nea 


(53) 


This field corresponds to the average of the function g carried by a fluid particle on site (i. j). The volu¬ 
metric coarse-graining is related to the areal one through 


(S)ij = 


h 8 ij 
hij 


(54) 


3.2.3 Definition of a velocity field on the grid 

Let us now define a large scale velocity field (or mean flow) on the uniform grid of the discrete model. 
We will introduce later a field accounting for small scale fluctuations of the velocity at each grid point. In 
the case of the actual shallow water model, for a given triplet of continuous fields h, q. II, the velocity 
field is computed by using Eq. (1201) . which involves two spatial differential operators, namely V 1 zW 1 and 
VZW 1 / 2 . Discrete approximations of these spatial operators are well defined on the uniform grid of the 
discrete model. Discrete approximations of Eq. (f20l> can therefore be used to define a velocity field on the 
uniform grid. Let us consider a coarse-grained vorticity field CO — h(q) — f and a coarse-grained divergent 
field jti = (ju) defined on the same uniform grid. Discrete approximations of the operators appearing in Eq. 
m can be written respectively as {V^A -1 [a>\} u = Lw Gf jkl .(O k i and {V^A -1 = Lu G^ kl .fi kh 

where the sum is performed over each grid site (k. /). In the remainder of this paper, we do not need 
the explicit expression of the kernels {G m } and (G' U } , which depend only on the domain geometry. Using 
these notations, we define the large scale velocity field as 

U mf,ij = {v^zr 1 \h(q) -f ]} + {vzi- 1 / 2 [(/x)]| , (55) 

where the index “mf” stands for “mean flow”. 

At this point one may wonder why the relevant coarse-grained fields used to define the large scale 
flow u m f should be h. (ju), ( q) (which, using Eq. ( l54l > . is equivalent to either the triplet li. hj. 1 . hq or to the 
triplet (h^ 1 ) , (id ), (q)). Our motivation for such a choice is twofold. First, we will see a posteriori that this 
allows to recover previous results derived in several limit cases (weak flow limit, quasi-geostrophic limit), 
and to obtain consistent results in the general case. Second, in physical space, the number of fluid particles 
at each grid site is given by the areal-coarse-grained height field h, according to the previous section. The 
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relevant macroscopic potential vorticity field or divergent field is then given by the volumetric coarse - 
graining (q) — hq/h and {/.() — hjl/h. In that respect, the phase space variable /? does not play the same 
role as q and u when considering macroscopic quantities in physical space, and this is why we do not 
consider the triplet ( h ), (u), (q) to describe the system at a macroscopic level. 


3.2.4 Definition of the micro canonical ensemble for the discrete model. 

Here we introduce a set of constraints associated with the discrete model in order to define the micro- 
canonical ensemble. These constraints are a discrete version of the Casimirs functional and the energy of 
the continuous shallow water model, defined in Eqs. (134b and ( 1321 ). respectively. Additional assumptions 
on the form of the energy will be required, and we will discuss the relevance of such assumptions. 

Note that we introduce here constraints for the discrete shallow water model, but we do not define what 
would be the dynamics of the discrete model. Indeed, it is not necessary to know the dynamics in order 
to compute the equilibrium state of the system. Only the knowledge of constraints provided by dynamical 
invariants is required. 

Let us consider a given microstates Xmicro — { (Jn ■ f i ) ■ c ln ■ 14 ■ /?„ } | > which belongs to the en¬ 

semble X m j cw of possible configurations of the discrete model defined in Eq. ( |481 ). By construction of the 
ensemble X mlcro , each element Xmicro satisfies the areal filling constraint given in Eq. (147b . 


The potential vorticity moments of the discrete model are defined as 


NM 


W - 0 ’ wl> 


72=1 


N‘ 


«r 


(56) 


The notation Y.ij means that the sum is performed over each grid point, with 1 < i < Nx, l<j< Ny. 
It is shown in Appendix |B] that those discrete potential vorticity moments tend to the potential vorticity 
moments of the continuous dynamics defined in Eq. (134b in the limit of large number of fluid particles AM. 


We also define the total energy of the discrete model as 


A = 


NM 




~ $cst ■ 


(57) 


where hj,n is the discrete topograph \Q. where M m f u is the mean flow defined in Eq. ( 155b . and where 
S' cst [hbij] is an unimportant functional of hi, chosen such that S — 0 at rest (i.e. when fi n = 0, u= 0, 
h n = 1 — hbi n j n for any n and M,-,- = M for any (/, j)): 


<%Cst ~ 



(58) 


A simple interpretation for this form of the total energy is that each fluid particle carries a kinetic energy 
associated with the mean flow u m f, as well as a kinetic energy associated with local fluctuations of the 
divergence field and finally a potential energy (the height of the center of mass of the fluid particle is 
hbi n j n T h n /2.). 

5 Here, the bottom topography is assumed sufficiently smooth to be considered constant over a grid site. A fluctuating 
topography would require further discussion. 
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It is argued in Appendix [B]that with only a few reasonable assumptions on the properties of the equilib¬ 
rium state, the energy of the discrete model defined in Eq. (157b would also be the energy of the equilibrium 
state of the actual shallow water model defined in Eq. (132b in the limit of large number of fluid particles 
NM. Note also that according to Eq. d57b . the vortical part of the velocity field does not contribute to local 
small scale kinetic energy, which is analogous to previous statistical mechanics results for non-divergent 
flow models such as two-dimensional Euler equations or quasi-geostrophic equations 12311331 . Qualita¬ 
tively, this is due to the fact that inverting the Laplacian operator smooth out local fluctuations of the 
relative vorticity CO — hq — 1 so that the streamfunction associated with the microscopic vorticity field is 
the same as the streamfunction associated with the coarse-grained relative vorticity field, see Appendix B 
for more details. 

The expression of the energy in Eq. ([571) involves a sum over the NM fluid particles. This sum can be 
recast into a sum over the N points of the grid, by using the definition of the areal coarse-graining in Eq. 
(150b and the definition of volumetric coarse-graining in Eq. (153b : 

^ ~ flY* ^2 ((^ My) + 2 + — l) + (j T ~ij ~ %) ) ^ • (59) 


This Eulerian representation of the energy allows to identify three different contributions. 

One first contribution to the total energy is given by the sum over each grid point {i. j) of the kinetic 
energy of the mean flow u,„/ (which is carried by M,j = Mhtj fluid particles), and of the potential energy 
of the areal coarse-grained height field hu : 


— ?v h ij U mf,ij + 8 (hij + hb,ij l) 


2 N 


(60) 


This contribution will be referred to as the total mean flow energy, or the energy of the large scale flow. 

A second contribution to the total energy is given by the sum over each grid points of the variance of the 
divergence levels fi carried by fluid particles at site (/, j) , times the number of fluid particles (M, ; = Mhjj ): 

^ ^ U l U (<M 2 ) y - My) • ( 61 ) 

This term can be interpreted as a subgrid-scale (or small scale) kinetic energy term due entirely to the 
divergent part of the velocity field (see Appendix iBb. 

The last contribution to the total energy is the sum over each grid point of the potential energy associ¬ 
ated with local fluctuations of the height field: 

(62) 

This term can be interpreted as a subgrid-scale (or small scale) potential energy term. 

One can finally check that total energy defined in Eq. d59b is the sum of the three contributions given 
in Eqs. (l60b . (l6~it and (l62b : 

& = $mf + <?5n + &5h- (63) 

It is well known that for the 2D Euler model or the quasi-geostrophic model, there is no contribution 
to the energy from the sub-grid fluctuations of the potential vorticity in the limit of vanishing grid size (see 
02211331 ). For the shallow water model, this is also the case. However the sub-grid fluctuations of height 
and divergence do contribute to the energy (see Eqs. ( 16Tb . ([62b and (l63b). A qualitative reason for this 
contribution of local height and divergence fluctuations to the total energy is that those two fields may 
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be decomposed on the basis of inertia-gravity waves, which are known to develop a small scale energy 
transfer II53I . and consequently to a loss of energy at subgrid-scales in the discretized model. 


Now that we have defined the configurations space in Eq. d48l) . the potential vorticity moments { \%micm} } k >o 
of the discrete model in Eq. d56l) and the energy S \%„,icro] of the discrete model in Eq. d59l ). we introduce 
the microcanonical ensemble as the restriction of the configurations space X m j cro to configurations with 
fixed values of energy E and Casimirs {Zf.} k>{) : 

{ Xmicro — {(Jn-Sn ) ■ tyn : } 1 <n<MN G Xmicro \ & [Xmicro] — ^G N S k \Xmicro \ — Z k | ■ (64) 

In our discrete model, we assume that the microcanonical measure over the ensemble of configurations 

Xmicro IS . 


, N,M 

a %,{Z*} t > o 


( Xmicro) 


<5 K [Xmicro] ~£) 

Q(E,{Z k } k > 0 ) 


I! S ([imicro] - Z k ) ]""[ 8 (1,, \% micro] - 1) 

k=0 ij 


NM 

X J | // .. 1 1// .; dty .; d u , ; d 1 n d./j, . 

n= 1 


(65) 


where Q (E. {Z^ }) is the phase space volume defined as 


N x N x 

Q(E,{Z k } k > 0 )= £ ••• £ 



h = 1 /*=!/!=! 


Jn= 1- 


r nm 


n h\dh n dq n djx n 


n= 1 


5(^[Xmicro] ~E) J] ^(^[Xmicro] ~Z k ) 
k=0 


X ^ ( by [Xmicro] 1) • (66) 

U 


The terms d/„ and d./ ( , are discrete measures with support on the grid coordinates: d/„(/) = Yjli fh 
d./„(/) = YjjL\ f j- The product ]~[;y is performed over the grid sites (i. j) with 1 < i < N x and l<j<N y . 

The constraint 8 (1 (/ [ Xmicro] I) corresponds to the area filling constraint defined in Eq. ( 15 11 1, which must 
be satisfied by each microstate Xmicro G X m i cro . 

Note that Eq. (f65T> is a discrete version of the formal microcanonical measure given in Eq. d40l > for the 
continuous case. The expectation of an observable ,vZ [ Xmicro ] ' n the microcanonical ensemble is 

(^)E,{Z k } k > o = d ^E’{Z k } k y o = j d V-E\z k } k > o (Xmicro ) && [Xmicro ] • ( 67 ) 

The problem is now to compute the macrostate entropy, which accounts for the logarithm of the number 
of microstates corresponding to the same macrostate. Although we do not need to consider this problem in 
the present work, it would be very interesting to have continuous or discrete approximation of the shallow 
water equation that have the invariant measure (165b . in an analogous way as what was achieved for the 2D 
Euler equations |[9lfT3l. 


3.3 Macrostates and their entropy 

The aim of this subsection is to compute the equilibrium state of the discrete model introduced in the 
previous subsection. The first step is to define the macrostates of the system, the next step is to compute 
the most probable macrostates, using large deviation theory, which yields also a concentration property 
asymptotically for large N and M (almost all the microstates correspond to the most probable macrostate). 
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3.3.1 Definition of the empirical density field 

Let us consider a microstate Xmicm — {(Ju-Ju) . c/,,. ./i,,} \< n <M!\ picked in the ensemble of configurations 

X ln icro defined in Eq. (148b . The empirical density field of this microstate is defined for each grid point (i. j) 
as 

Pij (o*, O q , Op) [Xmicro] = 8{h- o h ) 8 {q - o q ) 8 (p - Op).., (68) 

where the overline operator is the areal coarse-graining defined in d50l . This field contains all the statistical 
information of the system at the grid level. The constraint that each grid site (i. j) is covered by particles 
is given by Eq. (15 lb . which ensures the normalization: 


j dOhdOgda^ = 1 . 


(69) 


Let us consider a function g(h n .q n .p n ) depending on the height, potential vorticity and divergence 
carried by a fluid particle. Let us then consider the discrete microscopic field n —>• g n = g(h n ,q n .p n ) with 
1 <n< NM. Following Eqs. (BUI) and d68l . the corresponding coarse-grained field g i • is expressed solely 
in terms of the empirical density field pij: 

dOhdOqda^ g{Oh,Oq,G i i)pij (<7/,, <7 f y, <7^) . (70) 

If we consider for instance the function g(h.q.p) — q, then g n — q„ the microscopic potential vorticity 
field, and coarse-grained potential vorticity field q ■ is obtained by a direct application of Eq. (1701) : Z/ (/ = 

/ dO h dOqdOq GqPij. 

Importantly, the constraints { ;i?). [Xmicro] — -^Ar}*>o an d [Xmicro] — E defined in Eq. (1561 ) and ( 1591 ) de¬ 
pend only on the empirical density field pij (since they depend only on local areal coarse-grained moments 
of the different fields). The empirical density is therefore a relevant variable to fully characterize the system 
at a macroscopic level. 



3.3.2 Definition of the macrostates 

The macrostates are defined as the set of microscopic configurations leading to a given value p,j — pjj of 
the empirical density field: 


P — { Xiuirm £ X mlcro | V7. j pij \Xmicro\ — Pij } • (71) 

For the sake of simplicity, we make a small abuse of notation by denoting p both the macrostate defined 
in Eq. (1711 and the field p = {pij}. The values of the constraints are the same for all microstates within 
a given macrostate since they depend only on the local coarse-grained moments of the different fields, 
which remain unchanged for a prescribed empirical density field. The energy and the Casimirs, defined in 
Eqs. ( l57l) and (1561 ) respectively, have the same values for all the microstates within a single macrostate and 
will therefore be denoted by S’ [p] and { 33'^ |p] } k A) . 




20 


3.3.3 Macroscopic obsert’ables and empirical density 


Let us now consider an observable sX [Xmicm\ on the configuration space X m i cro defined in Eq. (148b such 
that their dependance on the microscopic configuration Xmicro occurs only through the empirical density 
field: 

■ l/ ^ [Xmicro] — &X [{ Pij [Xmicro] }] ■ (72) 

This is actually the case for any observable written as a sum over the fluid particles, i.e. for any observable 
appearing in Eq. ( 167b . It is therefore possible to change variables from Xmicro to the empirical density field 


values { Pij } 


1 <iJ<N x N y 


in Eq. (167b : 


W) 


NM 

E,{ Z k}k> 0 


n 




A(P) 

&{£,&}*> 0 ) 


8(£\p] 


E)f\8mp]-z k ) 

k=0 

x n 5 ( T 6 [P(./] - 1 ), (73) 

ij 


where 


N X N X Ny Ny nr NM 1 

fl(P)= £ " E L L / WdhndqndUn ^1 [^(P0'Xmicro] “ Pij)] (74) 


6 = 1 In=U 1 = 1 Jn = 1‘ 


- n= 1 


is the volume of a macrostate p defined in Eq. ( 171b in the configuration space X m j cw defined in Eq. (148b . and 
where £2 ( E,{Z k } k>0 ) is the total volume in phase space defined in Eq. d66t. The constraint 8 (1 ij [Pij] ~ I) 
is a normalization constraint for the {pij}, since 1, 7 |p (/ ] = f dc7/,df7 (/ d(7 [( pij. This normalization constraint 
also corresponds to the constraint that the fluid particles within a site must fit the available area, see Eq. 
(147b . The term f 6 /j [p, 7 ] means that the integral is performed over all the possible functions pij. The term 
8 {pij - pij) is a Dirac delta distribution on the functional space of the empirical density field values 
pij((T/i. <J q . (Jq ), with pij (oh,a q , <7 U ) [ Xmicro} the empirical density field defined in Eq. ([68b . 


3.3.4 Asymptotic behavior of the macrostates volume and derivation of its entropy 


Let us compute the asymptotic form of Q (p) defined in Eq. ( f74l i. by considering the limit M —>• where 
M is the average number of particles per grid site (/, j). For a given set of macrostates p defined in Eq. 

(ZD- the number of fluid particles per grid site is 

Mij — M J do h do q dGji Gi,Pij. (75) 

This is the only constraint on the particle positions in the grid for a microstate that belongs to the macrostates 
p. All the realizations of the particle positions that satisfies ( 1751) count equally in the ensemble of macrostates 
p. Through combinatorial calculation, the number of realizations of {(/„, J n )} that satisfy (1751 ) is (NM )! / IXy 37,/ 
Using Eq. (152b to express A7,y in terms of hij, we get 




(76) 


where 


Qjj — 


Mh: 


J Mhij J 


: h ,,,d// j-,,d(j tip,,, 8 I ^ . 8 (h m G/ t )8(q m G q ) 8 (p„, Gu j P// 


- m= 1 


m= 1 


(77) 
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is the number of possible configurations for a given set of M ,, = Mhjj fluid particles at site (/. j). 
The asymptotic behavior of the pre-factor in Eq. ( 1761 ) is computed through the Stirling formula: 
(NM)\ \ 


log 




( 78 ) 


The asymptotic behavior with M of Qjj defined in Eq. (1771) can be computed by using Sanov’s theorem. 

Before applying this theorem to our problem, let us consider the simpler case of K independent and 
identically distributed variables {Xk}\<k<K with comm °n probability density function F (x)- Those vari¬ 
able take values in a bounded interval of Sanov’s theorem describes the large deviation of the empirical 
density distribution 

1 K 

fK = ^Y. 5 (Xk~X), ( 79 ) 


k =1 


which can be considered as an actual probability distribution for the variable X- The probability distribution 
functional of this empirical density function is 


^\f\= [ f\F(Xk)dx k S(f-fK). 

7 k=\ 


( 80 ) 


Sanov’s theorem is a statement about the asymptotic behavior of the logarithm of &\f]- For a given 
function /(^), Sanov’s theorem states 


log mf\) K ^-K j dx/U)log(£M) 


( 81 ) 


if fdxf (x) — 1 ‘irid log (FX [/]) ~ — °o otherwise. An heuristic discussion of Sanov’s theorem is given in 

Ref 63. 

Combining Eqs. ( [79] ). (l80l > and ( ISTl ) and generalizing this result to K independent and identically dis¬ 
tributed L-tuple of variables \ { Xi.k } \ /. / f with common probability density function F ({Xi} kkz.)’ 

Sanov’s theorem is written in compact form as 


io § (/rp ({^i i</</.) n d %i,k & (/- ^ l n 5 j J K z~ k J n ^ ^ (j 

( 82 ) 

if/niid^i /(tel 1 < 1 < l ) — 1- 

Let us come back to the asymptotic behavior of Q,j defined in Eq. ( 1771 ). in the large M limit. Before 
applying Sanov’s theorem, it is needed to recast this equation into a form similar to the argument of the 
logarithm in the lhs of Eq. (fS2l i. Because of the 1 /h m term appearing in the delta Dirac function of Eq. 
dm one needs to perform first a change of variable from p (/ to 


— ° h „ 
K ij — y P‘ 
h u 


IJ 


( 83 ) 


6 


The fact that the variable %k have to be bounded is the reason why we set the cutoffs on the values of h n , q n and fi n . 






22 


The formal Jacobian ^ arising from this change of variable in the functional delta Dirac function depend 
only on hjj which does not depend on M and will therefore not matter for the asymptotic behavior of 
log(J2, ; ). In addition, the factor h 3 m appearing in front of the Lebesgue measure in Eq. ( 1771 must be 
divided by a normalization factor Z = f h 3 m dh m dq m dp m so that P{h mi q m ,Hm) — tf n /Z can be interpreted 
as a probability distribution function^ Then fl,j writes : 


Qq — 


■Mhi, 


Mh <i //,3 \ 

0 “ 7 ) dh m dq m dHm8 
m= 1 V Z / 



1 

Mhij 


MKij 

£ S ( h m 
m= 1 


O h ) & (^772 


A direct application of Sanov’s theorem to Eq. ([Ml then yields 



log hij 


j dOhdOqdOn Kij log 



(85) 


The normalization constraint on the distributions 7T (/ is already fulfilled by the definition of the coarse¬ 
grained height fields hq — lijj J d( 7 /,da f/ do'^ Kjj (see Eq. (1701) ). The inverse change of variable pjj — 
hi jTlj j j a h in Eq. ([85} yields 


log ( Q ‘j) ~ M 


j dGjjdGqdOjj Ghpij log 



( 86 ) 


Combining Eq. (1761 with Eqs. ([78} and (l 86 l . we obtain the asymptotic behavior of Q (p) with M: 


log (fl(p)) ~ MNy[p ], 


where 

y\ P ] 


is the macrostate entrop }0 



dG h dGqdGjj G h pij log 



(87) 


( 88 ) 


3.4 Continuous limit 


3.4.1 Expressions of the macrostate entropy, energy and potential vorticity moments 

Considering now the limit of an infinite number of grid site (N —> °°), the site coordinates (/, j) tend toward 
the continuous space coordinates x and the macrostate entropy y derived in Eq. ( 1881 ) becomes 


y\p] = - J dxdGhdGqdGji G h p (x, G h .Gq.Gj, ) log 


P (.A' 0/7 7 tJq■ Gq 

Gt 


(89) 


7 Here, Z depends on the cutoffs introduced in subsection l3.2l But we will see that it will vanish from the expression of 
the entropy in the end. 

8 Here we dropped the term log N coming from Eq. 1781 as it is constant that can be discarded by redefining Q E ,{z k } t>0 - 
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The empirical density has become a probability density function (pdf). For any function g : n^g n — 
g(hn,q n ,Hn)i its continuous coarse-grained field is now computed through 

do h do q do^ p (x, a h , o q , a^g (o h , o qi o M ), (90) 

The discrete mean flow defined in Eq. d55l) becomes (by construction): 

u mf \hM,hjI\ =V ± A- 1 (T^-f)+VA- 1 ' 2 (^J, (91) 

where all the coarse-grained fields have been expressed in terms of an areal coarse-graining, using Eq. 
d54). Similarly, the potential vorticity moments defined in Eq. (156b for the discrete model become 



VifceN 



The total energy defined in Eq. d59l i for the discrete model becomes 


(92) 


S 



*“«/ + 



+ g(h + h b - l) 2 


(93) 


We note that the total energy is a functional of li,h 2 ,hq,hp,lip 2 . Just as in the discrete case, this energy 
can be separated into a mean flow contribution (a large scale contribution including kinetic and potential 
energy), as well as a contribution from small scale kinetic energy due to local fluctuations of the divergent 
velocity field and a contribution from small scale potential energy due to local height fluctuations. The 
large scale (or mean flow) energy defined in Eq. (l60l > for the discrete model becomes 


S’mf \h,hq,hjl\ 


1 

2 


dx 


h u mf + g (h + h b — 1) 


(94) 


The small scale (or subgrid-scale) kinetic energy due to local fluctuations of the divergent part of the 
velocity field defined in Eq. (|6TT > becomes 




Kp,p 2 




(95) 


and the small scale (or sub-grid scale) potential energy due to local height fluctuations defined in Eq. ( 162b 
becomes 


&5h 




(96) 


One can check that the total energy in Eq. ( l93l ) is the sum of the three contributions given in Eq ( l94l ). ( l95l ) 
and ( 1961 ): 


$ — S'mf + S’sn + ^Sh- 


(97) 
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3.4.2 Microcanonical variational problem for the probability density field 


Let us come back to the average of a macroscopic observable .'// defined in Eq. (173k Using the asymptotic 
estimate for Q [p] given in Eq. (TsTb . the average of an observable /// defined in Eq. (Tub becometQ 


Win*# = I ®\P] p\ 


q,h,H 


£ NMy\p\ 

n(E,{z k }) 


8 [g [p] -E)Y\8 (& k [p] - Z k ) n 5 (1 [p] - 1). 


(98) 


k =0 


It comes to a Laplace-type integral where NM —> oo. Thus, the value of (///), m,n will be completely 

dominated by the contribution of the pdf p that maximizes the macrostate entropy defined in Eq. < f89l) 
while satisfying the normalization constraint 

T [P] = j dOhdOgdO^ p (X, Oft, o q , Op) = 1, (99) 

and the microcanonical constraints g\p]=E ,{3f k [p] = Z k } k> 0 , where the energy is defined in Eq. d93t 
and the potential vorticity moments are defined in Eq. ( 192b . This variational problem can be written in 
compact form as: 

max{^[p]| g\p] =E, VkeN S k [p]=Z k ,^€S T(x) [p] = 1} . (100) 

P 


The probability measure p induced by the empirical density field has a concentration property. In other 
words, the average of an observable depending only on macrostates is dominated by the most probable 
macrostates, which are solutions of the variational problem ( 1 100b . 

An interesting limit case for the entropy in Eq. d89b is worth mentioning in order to relate the variational 
problem in Eq. ( 1 1 00b with previous studies on the shallow water system. Let us assume that there is neither 
height variation nor divergent fluctuations, which would be the case if considering a quasi-geostrophic 
model or 2d incompressible Euler equations. Then the only height level and the only divergence level are 
Oh — \, — 0, respectively. Defining p q (x, of) — p (x, I. O q .0), Eq. ( f89l ) becomes up to an unimportant 

constant: 

[p ? (x, of} = - J dxdctq Pq log pq. (101) 

We recover in that case the macrostate entropy of the Miller-Robert-Sommeria theory 1221134b . 

Let us now assume that the height varies with position but that there is no local height fluctuations. 
Then at point x the only height level is a ), = h(x) — h(x). Defining p, /U (x, O q . ay,) = p (x,h(x),O q , Ou ), 
the macrostate entropy in Eq. (189b becomes up to an unimportant constant: 

y [Pqn (x,Oq,<Tn ), h (x)] = -J dxda^dau hpqn logp w . (102) 

This form of the entropy was proposed by Q, without microscopic justification. Interestingly, fTj obtained 
Eq. (1102b by assuming that the macrostate entropy can be written as 


y = 


j dxda^do,, h(x)s (p qtl (x, o q , o ^)), 


(103) 


and by noting that s(p) = — plogp is the only function that leads to equilibrium states that are stationary. 

9 Strictly speaking, the equal sign should be noted x which means that the logarithm of the terms on both sides are 
equivalent, see e.g. Ref. ED. 
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4 General properties of the equilibria and simplification of the theory in limiting cases. 


General properties of equilibrium states, solutions of the variational problem in Eq. (1 100b . are discussed in 
this section. Critical points of the variational problem are given in the first subsection ; they are computed 


in anncndixlcl This allows to obtain an equation for the large scale flow, and to show that equilibrium states 
of the shallow water model are positive temperature states. A weak height fluctuation limit is considered 
in a second subsection. It is found that the large scale flow and the small scale fluctuations are decoupled 


in this limit, and that there is equipartition between small scale potential energy and small scale kinetic 
energy. The quasi-geostrophic limit is investigated in a third subsection. 

4.1 Equilibrium states are stationary states with positive temperatures 

4.1.1 Properties of the critical points 

Critical points of the equilibrium variational problem defined in Eq. (11001) are solutions of the equation 



(104) 


k=0 J J 

where /3, {a*} t>0 and q (x) are the Lagrange multipliers associated with the conservation of the energy. 


the potential vorticity moments, and with the normalization constraint, respectively. The computation of 
the critical points is performed in Appendix Icl 

A first key result of Appendixes that solutions of Eq. (11041 ) factorize as 


p=p h (x, a h ) p q (x, a q ) p^ (x, Off). 

A second result of Appendix O is that the pdf of the divergence field is a Gaussian 


(105) 



(106) 


Recalling that j3 = dS/dE is the inverse temperature, we see that the variance of the pdf of the divergence 
field is given by the temperature of the flow: 


— ~jj~ = i. 


(107) 


This equation has important physical consequences. First, fluctuations of the divergent field do not vary 
with space. Second, the temperature of the equilibrium state is necessary positive. This contrasts with 
equilibrium states of two-dimensional Euler flow, which can be characterized by negative temperature 


states Q3II231 . This was realized first by Onsager for the point vortex model lf28l . In the context of 2D 


Euler equations in doubly periodic domains, the equilibrium states are always characterized by negative 
temperature, just as in the case discussed by Onsager (this result is shown in @). However, the existence 
of large scale flow structures characterized by positive temperatures are possible at low energy in the 
presence of lateral boundaries and non-zero circulation, and/or in the presence of bottom topography in 
the framework of QG model. Such positive temperature states are known and documented in the literature 
(see e.g. the original papers by li23l , j34l and fl8l or j4j and references therein). 

In addition, it follows from Eqs. ( l95l ). (11051 ) and ( 11071) that the temperature is directly related to small 
scale kinetic energy due to the divergent velocity field: 



(108) 
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A third result of Appendix ICl is the expression of the pdf of the height field: 


Pi i (*, o*) = Q-^y a /T ex P ( -P'kGh - 


$h (*) 

Oh 


G h (x) = j d a h of exp ^-j3| o A - 


& (*) 

Oh 


where <Sj/, (x) a function related to h (x) through /t = / do/, 0/,P/,. 

A fourth result of Appendix[C]i s the expression of the pdf of the potential vorticity field: 


(109) 


p q (x,o q ) = 


G q {pW mf ) 


exp pVmfCJq - £ VkOq , G 9 ([3%,f ) = / da 9 exp P%,fO q - £ a k a k q 


k= 1 


*=1 


( 110 ) 

where W m f \Ti. q. pj is the mass transport streamfunction of the mean flow defined through an Helmholtz 
decomposition of ha m f. 

7m mf = V L % lf + V<P lnf . (Ill) 

According to Eq. (IllOb . the coarse-grained potential vorticity field is a function of the mass transport 
streamfunction: 


r— 9 


G'(pW mf ) 


’ Gq(pw mf y (112) 

A fifth result of Appendixes that the mass transport potential of the mean flow defined in Eq. ( 11111) 


vanishes: 

and the velocity field can now be written 


^mf — 0, 


(113) 


u mf = =V ± W mf . (114) 

h 

A sixth result of Appendix [Cl concerns the Bernoulli function of the mean flow defined as 

B mf = ] -u 2 mf + g (Ti + h h — I). (115) 

According to Eq. (1218b . combining Eq. ( 12161 ) and Eq. ( 12171 ) of AppendixICl allows to express B m f in terms 
of [5, Gq, Gh and a constant Aq that can be computed in principle using the conservation of the total mass 

*> = 1 : 

B m f — P 1 log(GgG^) +g/z+Ao. (116) 


4.1.2 Equation for the large scale flow 

Let us now establish the equations allowing to compute the large scale flow. A first equation is obtained 
by injecting the expression for the mean flow u m f given in Eq. (II 141) into the expression of the Bernoulli 
function defined in Eq. ( II 15b . and by combining Eq. (II 15b with Eq. (II 16b : 

+ghb = o-log(G 9 G/,) +Ai, (117) 

1 h P 

with Ai = + g. A second equation is obtained by taking first the curl of Eq. (|9H . which yields hq — f — 
V x u m /. Then, replacing uby its expression given in Eq. (II 14b . remembering that the potential vorticity 
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field and the height field of the critical points of the variational problem are decorrelated (qh = hq), and 
replacing q by its expression given in Eq. (II121) yields 


T & q (Wnf) 

(P %,f) 


/ = V 



(118) 


The closed system of partial differential equations ( 1117b and (1118b must be solved for h and W m f for a 
given value of [1 , G q , G/„ and A\. The set of parameters [5. { M i must in fine be expressed in terms 

of the constraints of the problem given by the energy E, and the potential vorticity moments {Z^} jt>0 . This 
may require a numerical resolution in the general case. 


We have seen in subsection 12.21 that a flow described by (h. u ) (or equivalently by (h.q.ii), or by 
(/;, <f>. W)) is a stationary state of the shallow water model if and only if J(q. W) — 0 and <P ~0. 

According to Eqs. (II 12b and (II 13b . the large scale flow (li,xx m f) is a stationary state of the shallow 
water model. 


4.1.3 Comparison with previous results 

The set of equations (II 17b and (II 18b describing the large scale flow is similar to the one obtained by 
Weichman and Petrich ll56l (through a Kac-Hubbard-Stratonovich transformation) and by Chavanis and 
Sommeria J7] (through a phenomenological generalization of the Miller-Robert-Sommeria theory), ex¬ 
cepted that the rhs of Eq. ( II 17b contains an additional term jr log (G/;) + A \ which does not appear in these 
previous works. The reason for the presence of this additional term is that we have taken into account the 
presence of small scale fluctuations of height and velocity, which were neglected in Refs. 11561171 . We will 
see in subsection 14.31 that this additional term becomes negligible with respect to the others in the quasi- 
geostrophic limit, in which case we recover exactly the set of equations for the large scale flow obtained 
in Refs. f56ll7l . 

In addition, we have shown in subsection ld. 1.1 I that only positive temperature states are allowed, which 
shows that only one subclass of the states described by Eqs. ( II 17b and (II 18b are actual equilibrium states. 


4.2 Equipartition and decoupling in the limit of weak local height fluctuations 

We consider in this subsection the limit of weak local height fluctuations. This step makes possible com¬ 
putation of explicit solutions of the variational problem ( 1 100b . Meanwhile, it allows to explore the conse¬ 
quence of the presence of these small scale fluctuations on the structure of the large scale flow. By "limit 
of weak local height fluctuations”, we mean 

Vxe (h 2 -h ) <C h. (119) 

As already argued in B561 . this limit of weak local fluctuations is physically relevant, since the presence 
of shocks in the actual dynamics tends to dissipate small scale fluctuations of height or kinetic energy. This 
will be further discussed in subsection [572] 
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4.2.1 The height distribution 


Assuming in addition that height levels <7/, such that (7/, 


h\ > 



do not contribute signifi¬ 


cantly to the pdf pi, defined in Eq. ( 11091 ). we perform an asymptotic development with (a*/A — l) <C 1, 
and obtain at lowest order in this small parameter a Gaussian shape for the pdf: 


Pli (X, <7/,) 



( 120 ) 


Similarly, the term G* defined in Eq. d 109b can be computed explicitly in this limit. One gets at lowest 
order 

G h = Ti (x) exp {-gfih (x)). (121) 

Using Eq. (1 1 20b . the weak local height fluctuation limit given in Eq. ( 11191 ) can be interpreted as a low 
temperature limit gh 1//3. Injecting Eq. (11211) in the previously established relation in Eq. (II 161 ) yields 


B m f = j log G q + ^ log h + A 2 , (122) 

where A 2 is a free parameter determined by the conservation of the total volume. In the rhs of Eq. (1 122b . 
there remains an additional term log/; which is not present in the large scale flow equation obtained 
by If56f 71, but considering the weak flow limit allowed to obtain an explicit expression for this additional 
term. We will see in subsection l4.3l that this term becomes negligible with respect to jr lcgG^ in the quasi- 
geostrophic limit. 


4.2.2 Equipartition and the limit of weak small scale energy 

Injecting Eq. (11 20b in the expression of the small scale potential energy defined in Eq. (1961 ) yields 

= (123) 

Comparing this result with Eq. ( 11081 ) shows equipartition of the small scale energy between the potential 
energy and the kinetic energy. The total energy due to small scale fluctuations is 

Efiuct = <?5h + &8n — -p- (124) 

This equipartition result of the small scale energy was already obtained by Warn when computing 
equilibrium states of the Galerking truncated dynamics in a weak flow limit l53l . More precisely. Warn 
decomposed the dynamics into vortical modes and inertia-gravity modes (which are defined as the eigen- 
modes of the linearized dynamics), and concluded that the energy of the equilibrium state should be 
equipartitioned among inertia-gravity modes in the limit of infinite wavenumber cut-off. This equipar¬ 
tition of energy among inertia-gravity modes would lead to equipartition between potential and kinetic 
energy at small scales, just as in our case. Here we have recovered this result with a different approach and 
we have generalized it beyond the weak flow limit. 
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Finally, we remark that the equipartition result of Eq. (11241) shows that the low temperature limit 
gh 3> 1 //3 corresponds to a limit of weak small scale energy due to local fluctuation of the height field and 
of the divergent field: 


^ Efluct • 


(125) 


4.2.3 Decoupling between the large scale flow and the fluctuations 


Still by considering the weak local height fluctuation limit, let us assume that Efi uct defined in Eq. (11241) 
is given, which means that the temperature is given. According to Eq. (II 181 ) and Eq. (1117b . knowing the 
coarse-grained height field h and the pdf of potential vorticity levels p 9 (x, o q ) (which allows to compute 
G^) is sufficient to determine the mass transport streamfunction f ' m f, and hence the velocity u m f by using 
Eq. (11 141) . Then the mean-flow energy defined in Eq. ( l94l ) does not depend on p^: 


&mf Pb Pq] — 


dx 


hu 2 mf [h,q] +g(h + h b -iy 


( 126 ) 


The total energy defined in Eq. (1971) is therefore the sum of the mean-field energy associated with a large 
scale flow that depends only on h and p q , and of the energy of small scale fluctuations associated with 
height and divergence fluctuations: 

& — &mf \fli Pq] T Epi UCt • (127) 

In addition, injecting Eq. (1 105 b in Eq. ( 1921 ) gives the expression of the potential vorticity moments as a 
functional of h and p q only: 


VleN 3T k \h, p q \ = j dxda, hp q a q . (128) 

Let us now consider the macrostate entropy functional defined in Eq. (189b . Injecting Eqs. (1106b and 
(1120b in Eq. (1105b . and performing the asymptotic expansion of the integrand of the mean-field entropy 
with (<J b /h — l) <C 1 leads at lowest order to (up to a irrelevant constant) 

& = #rnf \h,p q } + Sfluct (■ E fluc ), (129) 


y m f pi, Pq] = - J dxda 9 hp q log ^ 2 , (130) 

Sfluct ( Efiuct ) = log (Ef l UCt ). (131) 

Since p q and h are two fields allowing to compute the large scale flow of energy £ m f, and since the 
fluctuations of the potential vorticity field do not contribute to the small scale energy Ef[ uct , the entropy 
S^mf [h. p q ] will be referred to as the macrostate entropy of the large scale flow. 

The second contribution to the total entropy in Eq. ( 1129b depends only on the energy of the small 
scale fluctuations Efi uct . Since this energy is solely due to the local variance of the height field and of the 
divergent field, it will be referred to as the entropy of the small scale fluctuations. Note that in that case 
the height field is involved both in the large scale flow through its local mean value, and in the small scale 
fluctuations through its local variance. 

The decoupling of the energy and the macrostate entropy functional into a part that depends only on 
p q ,h and another part that depends only on small scale height and divergent fluctuations with energy Efj uct 
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has both a useful practical consequence and an interesting physical interpretation. The variational problem 
in Eq. (1 100b can now be recast into two simpler variational problems: 

S [E , {Z k } k>Q ) = max { S m f [E — Efj ucl , + Sfj ucl (Efi uc t) } j (132) 

Efluct 

Smf(E mf ,{Z k } k > 0 ) = max _{^ mf \h,p q ] \ S mf [h,p q ,] = E mf , VteN 2E k [p q ,h\ = Z k } , (133) 

where ■9'mf, £ m f and are the functional defined in Eqs. (1130b . ( 1126b . and (1128b . respectively. The 
variational problem in Eq. (1132b describes two subsystems in thermal contact. In order to compute the 
equilibrium state, one can then compute independently the equilibrium state of each subsystem, and then 
equating their temperature in order to find the global equilibrium state. This classical argument follows 
directly from the maximization of (1132b . If the two subsystems can not have the same temperature (for 
instance when the temperature of both subsystem have a different sign), all the energy is stored in the 
subsystem with positive temperature in order to maximize the global entropy. 

In the present case, a first subsystem is given by the large scale flow of energy E m f, which involves 
the field h and the potential vorticity field described by the pdf of vorticity levels p q . The equilibrium state 
of this subsystem is obtained by solving the variational problem in Eq. (1133b . This variational problem 
corresponds to the one introduced by Q except that large scale flow energy is not the total energy but 
only the available energy when the energy of the fluctuations has been removed. The entropy of the large 
scale flow (1130b in the variational problem (1133b is closely related to the entropy introduced in Q (up to a 
functional of h). The potential vorticity moment constraints apply only to this subsystem. These additional 
constraints are essential since they allow for the possible existence of a large scale flow, see e.g. SI- 

The second subsystem is given by the local small scale height fluctuations and the local small scale 
divergent fluctuations with total energy Efi uct , and with the entropy given in Eq. (1131b . The inverse temper¬ 
ature of this subsystem fi — dS fiuct /^^-fiucu which, using Eq. (1131b . yields a temperature of (5 1 = Efi uct . 

In practice, it is easier to compute directly equilibrium states of the large scale flow subsystem in 
the canonical ensemble, where the energy constraint is relaxed. In that case the equilibrium states of 
this subsystem depend on the temperature Efi uct and on the dynamical invariants {Z k } k>{ y One just need 
to check a posteriori that this ensemble is equivalent to the microcanonical one by verifying that each 
admissible energy E m f is reached when varying Efi uct form 0 to -\-°° , for a given set of potential vorticity 
moments {Z k } k> Q . In order to find the actual equilibrium state associated with the total energy E, one 
then needs to solve the equation 

E — E m f ( Efi uct ,{Z k } £>o) + Efi UCt • (134) 

To conclude, it is now possible to study independently the large scale flow subsystem, to consider if 
necessary any approximation on this flow, such as the quasi-geostrophic limit or the Euler 2D limit, and 
finally to couple this subsystem with the small scale height and divergence fluctuations subsystem in order 
to select the actual equilibrium state. If one linearize the large scale flow entropy entropy (1130b and the 
large scale flow energy (1126b . this picture of two subsystems in thermal contact gives a justification to the 
variational problem introduced by 1451 and extended by ll25l . 14511251 suggest that for a given frozen in 
space potential vorticity, the dynamics should relax through geostrophic adjustment to a state minimizing 
the total energy. This result is recovered from the coupled variational problem (1132b . 

4.2.4 Either a non-zero circulation or a non-zero bottom topography is required to sustain a large scale 
flow at equilibrium 

It is shown in appendix [D] that when circulation is zero (i.e. when Z\ — f) and when topography is zero 
(h b — 0), a state with a large scale flow at rest (u m f = 0) is a maximizer of the large scale flow macrostate 
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entropy among all the possible energies: 

S(0, {Z k } k > 0 ) = max {S mf (E mf , {Z k } k > 0 )} when Z\ — f and h b = 0. (135) 

E, n f 

According to Eq. (11311) . the fluctuation entropy Sfi uc t increases with the fluctuation energy Efi uct . The total 
macroscopic entropy S m f (E — Efi uct ,{Z k } k>0 ) +Sfi uct (Efiuct) is therefore maximal when all the energy 
is transferred into fluctuations (Efi uct = E). This generalizes to a wider range of flow parameters and to 
a wider set of flow geometries a result previously obtained by Warn in a weak flow limit for a doubly 
periodic domain without bottom topography l53l . In addition, we will see in the next section that when 
there is a non-zero bottom topography and rotation, a large scale flow can be sustained at equilibrium. 


4.3 The quasi-geostrophic limit 

We show in this subsection that the variational problem of the Miller Robert Sommeria theory is recovered 
from Eq. (11331 ) for the large scale flow subsystem when considering the quasi-geostrophic limit, which 
applies to strongly rotating and strongly stratified flows. 


4.3.1 Geostrophic balance 

Let e}J} be the typical velocity of the large-scale flow and let L = \J\EZ\ be the typical horizontal scale 
of the domain where the flow takes place. We introduce the Rossby number and the Rossby radius of 
deformation respectively defined as 


Ro = 



f 


( 136 ) 


Here / is the Coriolis parameter, H — 1 is the mean depth and g the gravity, see subsection 12.11 If f* o 
we can always rescale time unit so that f — 1, and we make this choice in the following. It is also assumed 
that the aspect ratio of the domain where the flow takes place is of order one, so that L = 1 since we chose 
length unit so that \2>\ = 1. The quasi-geostrophic limit corresponds to small Rossby number, and to a 
Rossby radius of deformation that is not significantly larger than the domain length scale: 


1 , R~ l = &{ ( 1 ). 


( 137 ) 


By construction, the mean flow is of the order of the Rossby number: |u„,/| ~ Ro. The coarse-grained 
interface height is given by 


rf = h-l+h b . ( 138 ) 

Let us assume that the spatial variations in fluid depth are small compared to the total depth El = 1 , 
with the scaling 77 ~ R -Ro. At lowest order in Ro , the mean flow Bernoulli potential defined in Eq. (11151) 
becomes B m f — R 2 r\. Remembering that we consider in addition to the quasi-geostrophic limit a weak 
fluctuation limit given by Eq. (1 125b . which can be expressed as R 2 S> J3 -1 , Eq. ( 11221 ) yields at lowest order 

B,nf — R 2 77 = j 3 “ 1 log G q + cst. 


(139) 
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This equation implies d6 m ^/d*f^„/- = q, consistently with what we expected for a stationary flow in the 
absence of small scale fluctuations, see subsection 12.21 Taking the curl of Eq. ( 1 1 39b and collecting the 
lowest order terms yields geostrophic balanccT 7 ] 

R 2 \ /L rf = u m/ . (140) 

Eq. (1140b also shows that the scaling hypothesis for Tf is self-consistent. 

It is remarkable that equilibrium statistical mechanics predicts the emergence of geostrophic balance. 
We stress that those results, which have been obtained through the introduction of a semi-Lagrangian 
discrete model, are valid whatever the amplitude of bottom topography variations (i.e. beyond the usual 
approximation h/, ~ Ro required to derive the quasi-geostrophic dynamics). By contrast, in the framework 
of a Eulerian discrete model, one would find that the large scale flow is not at geostrophic equilibrium 
unless hb ~ Ro or h/, <C Ro, see Appendix E. 

4.3.2 Quasi-geostrophic dynamics 

The geostrophic balance is not a dynamical equation. When hb ~ Ro, the dynamics if given by the quasi- 
geostrophic equations. At lowest order in Ro, we get hu m f = u m f, and 

Vmf = i ^ l nj ■ Qtnf — ^ rn J — ( 141 ) 

where W m f and \ff m f are the transport streamfunction and the streamfunction obtained through the Helmholtz 
decomposition of hu m f and u m /, respectively. In that case, the relative vorticity is 

65 =A\jf m f. (142) 

The geostrophic balance <1 140b is equivalent to Tf — R 2 Ymf + C. The value of Xf/mf at the domain 
boundary can always be chosen such that the constant term vanishes, which yields 



Mass conservation given in Eq. (O leads then to the following constraint on the streamfunction: 

/ dx \f/ mf = 0. (144) 

Given a potential vorticity level o q , a change of variable can be performed by introducing quasi- 
geostrophic potential vorticity levelJH 

O g =(e q -\). (145) 

The pdf of quasi-geostrophic levels is 

Pg (x : Gg) = Pq (x, 1 + Og). (146) 

10 For the shallow water model, the fluid is at hydrostatic balance. Thus the pressure in the fluid is P(x,y,z,t ) =P 0 + pg(H + 
r](x,y,t) — z). Then the pressure horizontal gradient is simply proportional to the interface height horizontal gradient. Hence, 
the geostrophic balance simply writes i? 2 V ± rj‘ = u 

11 This change of variable is a guess guided by the fact that the QG potential vorticity is usually obtained by expanding 
the SW potential vorticity in the limit of small height variations minus a constant and unimportant term. Here we start by 
removing the unimportant constant (1 in our unit system) from the potential vorticity levels, and then perform the small scale 
expansion in height. 
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The local quasi-geostrophic potential vorticity moments are defined as 

WkeN,q k g = jdo g o g p g . (147) 

At lowest order in Ro, the coarse-grained quasi-geostrophic potential vorticity obtained by considering 
k= 1 in Eq. ( 1 147b becomes 


q g = ( 0 -T)+h b . (148) 

which, using Eqs. (1 142b and (1143b . yields 

q g =AYmf-^ + h b . (149) 

In the quasi-geostrophic limit, the large scale flow is fully described by the streamfunction I j/ m f, which can 
be obtained by inverting the coarse-grained potential vorticity field defined in Eq. (1149b . 


4.3.3 Quasi-geostrophic constraints and variational problem 

At lowest order in Ro, and after an integration by part, the expression of the mean-flow energy in Eq. (1126b 
is equal to the quasi-geostrophic energy: 


$mf,g 



(yy mf ) 2 + 


R 2 ’ 


(150) 


where t j/ m f can be expressed in terms of q„ through Eq. (1149b . Similarly, the conservation of the potential 
vorticity moments defined in Eq. (1128b implies the conservation of quasi-geostrophic potential vorticity 
moments 


VkeN, 3T g _ k = I dx q k g , (151) 

and the macrostate entropy of the large scale flow defined in Eq. (1130b writes now 

ZZmf,g Pg I d\dO\, p g log Pg. (152) 

Using Eq. (1150b . (1151b and (1152b . the variational problem in Eq. (1133b is now recast into a simpler varia¬ 
tional problem on the pdf p g : 

S m f, g ( E mf , {Z k } k> j) = max {y m f, g [p g ] \ & m f, g [pg] = E mf , VkeN [p ? ] = Z k } . (153) 

PgJPg =1 


The entropy the energy S m f g and the potential vorticity moments { :Z g k }, > : are defined in Eqs. 

(1152b . ( 1150b and (1151b . The variational problem defined in Eq. (1153b is the variational problem of the 
Miller-Robert-Sommeria statistical mechanics II22II33I . 
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4.3.4 Maximum energy states and consistency of the quasi-geostrophic approximation. 

We have shown that in the weak height fluctuation limit and the quasi-geostrophic limit, computation of 
the large scale flow associated with the equilibrium state amounts to the computation of the solution of the 
variational problem in Eq. (11531) . with the restriction that the temperature is positive (due to the coupling 
with small scale fluctuations of height and divergence, as discussed in subsection l4.1.1l ). Here we discuss 
the solutions of this variational problem, which are energy maxima for a given set of potential vorticity 
moments {3fg,k} k>l - Although the initial weak fluctuation or quasi-geostrophic limit may not be fulfilled 
for such states, they can always be computed, and it provides an upper bound for the energy of the large 
scale flow obtained in those limits. 

It is known that for positive temperatures (j3” 1 = Efi uct > 0), the equilibrium entropy defined in Eq. 
(11531 ) is concave 01, and the energy increases when /) decreases. This means that the state with a maxi¬ 
mum energy is reached when /3 —^ 0, see e.g. (331 . 

Injecting /) = 0 in the expression of p q in Eq. (11101) leads to a uniform mean potential vorticity field. 
According to Eq. (11461 ) and (11471) . this implies that the quasi-geostrophic potential vorticity field is also 
uniform: q g = Z g \ where Z g \ — f dxq g is the circulation. Since 7j„ is a constant, this state is referred to 
as the “mixed” state. We get 

Ayr mix -^-=Z gi i-h b . (154) 

Let us call E mlx the energy of the mixed state: 

E mix = max { S mf \\/k G N 2? gtk = . (155) 

PgJPs =1 

Using Eq. (1 1 5()t . it can formally be written 

Emix = ~\j ( Z sd “ h b) - Jp ) ( Z ?4 “ h b ) ■ (156) 

For a given domain geometry, a given circulation Z g \ and a given bottom topography field hi,, a non 
trivial large scale flow can be observed whenever E m i x > 0. We see from Eq. (11561 ) and (11551) that the con¬ 
dition for a large scale flow to exist is that either the circulation Z gi \ is non zero or the bottom topography 
hi, is non-zero. If both Z g \ = 0 and hi, = 0, then the energy of the large scale flow vanishes (E m f — 0), 
and all the energy is lost in small scale fluctuations (E — Efi uct ), consistently with the results of subsec¬ 
tion 14.2.41 In this case, coupling thermally a large scale flow with fluctuations leads to a state with all the 
energy lost in fluctuations. Note that in the case Z g \ — 0 R ~ 1, E m \ x can also be interpreted as a norm of 
the topography field hb- 

Since E m f < E m j x and since E m j x depends only on the problem parameters (namely the circulation, the 
Rossby radius and the bottom topography), a sufficient condition to have Ro <C 1 is 

4i 2 « 1- d57) 

If this condition is fulfilled, the quasi-geostrophic assumption is self-consistent (as well as for the scaling 
hb Ro in the case Z g \ — 0). 
























35 


5 Explicit computation of phase diagrams and discussion 

The aim of this section is to apply the results of the previous section to the actual computation of equilibria 
and their energy partition. In order to solve analytically the variational problem of the statistical mechanics 
theory, we focus on a subclass of equilibria referred to as the energy-enstrophy equilibrium states. This 
allows to build phase diagrams in a two parameter space, and to discuss the energy partition between a 
large scale flow and small scale fluctuations when these parameters are varied. We finally discuss the role 
of shocks that occur in the actual shallow water dynamics, and present a geophysical application to the 
Zapiola anticyclone. 

5.1 Energy-enstrophy equilibria for the quasi-geostrophic model. 

We consider the variational problem 

Sg,mf {E,Zj) = max {dZ 7 m f ,g [Pg] | <*mf,g [Pg] — E m f, [Pg] — ^2 , -2gi [Pg] — 0 } > (158) 

PgJPg= l 

where the functionals y m f,g, 3‘mf.g, »Zg 1.2 are defined in Eqs. ( 11521) . (1 1 50b . and ( 1151b . respectively. 

The peculiarity of this variational problem is that only two potential vorticity moments (the circulation 
and the enstrophy) have been retained as a constraint, in addition to the energy. Such energy-enstrophy 
equilibria are a subclass of statistical equilibria solutions of the more general variational problem given in 
Eq. ( 1153b . see e.g. |[Tll26l. For a given global distribution of potential vorticity, several limit cases on the 
energy allow to simplify the computation of the solutions of 

the variational problem in Eq. (1153b into the computation of the simpler variational problem in Eq. 
( 1158b . For instance, assuming that bottom topography is non zero, and that the global potential distribution 
is such that the mixed state q = cst exists, the solutions of (1158b are the solutions of the more general 
variational problem Eq. ( 1153b when E —> E m j x . 

The set of all potential vorticity fields q corresponding to solutions of the variational problem in Eq. 
(1158b have been previously described by lt6] l49ll5Tll26l in the case of a bounded geometry, and phase di¬ 
agrams were obtained with energy E and circulation Z\ as external parameters. The role of enstrophy Z 2 
was not discussed. The main reason is that for a given large scale flow characterized by E\ and Z\, chang¬ 
ing Z 2 would only imply changes in the small scale fluctuations of potential vorticity levels, see e.g. 1261 . 
and such small scale fluctuations do not contribute to the total energy. In the present case, Z 2 will play an 
important role in determining the energy partition between the large scale (vortical) flow and small scale 
fluctuations due to the height and divergent velocity field. For the sake of simplicity we consider vanishing 
circulations Z\ — 0. 

The problem 1158b is solved in Appendix [F] for positive temperature states, and we present here the 
main results. A typical phase diagram is shown on Fig. [2] This figure is obtained by assuming that the 
bottom topography is proportional to the first Laplacian eigenmode (see Fig. Oc) but it is explained in 
Appendix |F] that this phase diagram is generic to any bottom topography^ 

There are two important quantities related to the height field: the maximum allowed energy E mix de¬ 
fined in Eq. 1156b . and the available potential enstrophy 

Z b = Jdxh 2 b , (159) 

12 For such a bottom topography, the topography, the stream function and the potential vorticity field are all proportional 
to each other for any initial condition for the enstrophy Zi and the energy E. That is why we do not show plots of the fl ows 
for different point of the phase diagram in Fig. [2] We rather choose to consider the case of the Zapiola drift in subsection l5.3l 
to see the effect of different value for the initial energy. 
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which is the maximum reachable value for the macroscopic enstrophy f dx , see Appendix [F] Both Z/, 
and E m i x are a norm for the height field. 

The phase diagram of the quasi-geostrophic energy-enstrophy ensemble restricted to positive temper¬ 
ature states is presented in Fig. [2}a. As explained in subsection 14.31 the energy of the large scale flow 
E m f can not exceed the value E mlx defined in Eq. (11551) . due the restriction of positive temperature states. 
Depending on the sign of Z 2 — Z/,, where the potential enstrophy Z/, is defined in Eq. (1 159b . the system 
behaves differently: 

- When Z 2 > Zj, the minimum admissible large scale flow energy is zero. 

- When Z 2 < Z b , there exists a minimum reachable large scale flow energy E m i„ (Zq) below which there 

is no equilibria. The curve E m , n increases from 0 to E m i x when Z 2 decreases from Z/, to 0. 

The thick black line in Fig. [2}a delimits the domain of existence for the equilibria. The thin dashed black 
curves represent the isotherm, i.e. the points of the diagram with the same value of Efi uct . Note that there 
is no bifurcation in this phase diagram; to each point (£)„/-, Z 2 ) corresponds a single equilibrium state, 
whose expression is given explicitly in Appendix[F] The structure of an equilibrium large scale flow above 
a topographic bump at low and high energy are presented in the last subsection. 

The phase diagram in Fig. fallows to discuss the energy partition between small scale fluctuations and 
large scale flow when the quasi-geostrophic flow is coupled to small scale fluctuations of the height field 
and divergence field (through the shallow water dynamics). The motivation behind the works of 8531 and 
1561 was the prediction of energy partition between large scale flow and small scale fluctuations. Here we 
provide for the first time an explicit expression for such energy partition. 

We have seen that in the weak height fluctuation limit, the temperature of the large scale flow at 
equilibrium is given by the energy of the small scale fluctuations of height and divergence fields. The 
parameters are now the total energy E and enstrophy Z 2 . The mean flow energy E m f and the fluctuation 
energy E fi uct are found by solving 


E — E m f (Ef luctl z 2 ) +Efi UCtl ( 160 ) 

which is easily done graphically using the diagram of Fig. [2]- a, and performed numerically in practice, 
see Appendix |F] for more details. The ratio of the large scale energy E m f over the total initial energy E 
as a function of the total initial energy E is shown on Fig.U-b for different values of initial enstrophy Zi. 
According to computations performed in Appendix |FJ four cases for the energy partition in the low energy 
limit are distinguished depending on the sign of Z 2 — Z/, and the scaling of Zi — Z\, with E : 

- When Z 2 < Zb, the ratio E m r/E tends to 1 when E tends to the minimal admissible energy (Z 2 ). 

- When Z 2 Z\y, the ratio E^q JE tends to 0 when E tends to zero. 

- When Z 2 > Zb with Z 2 — Z/, ~ E a with a > 1/2, the ratio E m f /E tends to 1 /5 when E tends to zero. 

- When Z 2 > Zb with Z 2 — Z/, ~ E 1 / 2 , the ratio E m t/E tends to a finite value (depending on the propor¬ 
tionality coefficient between Z 2 — Z/, and E i 2 and the bottom topography). 

- When Z 2 > Z/, with Z 2 - Z/, ~ E a with a < 1/2, the ratio E m f/E tends to 0. 

Whenever Z 2 > Zj, we found Efi uct ~ E. We see on Fig. [2]-b that the ratio E m f/E converges to one when 
E —» E m in (Z 2 < Zb), that it converges to zero for Zi > Zb, and that it converges to 1 /5 for Z/, = Z 2 . Note 
also that E m f is bounded by E mlx such that E m f/E tends to zero when E tends to infinity. 
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Fig. 2: a) Phase diagram of the energy-enstrophy ensemble in the plane (Zt. E m f). Z* is the maximum 
reachable value for the macroscopic enstrophy defined in Eq. (11591) . E m i X is the maximum reachable energy 
for the mean-flow defined in Eq. (11551 ). The dashed lines corresponds to isotherms, on which the energy 
Efine of the small scale fluctuations is constant. The thick black line on the bottom left corner is a boundary 
below which no equilibria exist, b) Ratio between the energy of the large scale flow E m f over the total 
energy E m f/E as a function of the total energy E = E m f + Efi uct . The different curves correspond to 
different values of the initial microscopic enstrophy Z 2 represented by horizontal marked lines on panel a. 
The x-axis is on a logarithmic scale, c) Colormap plot of the bottom topography /?/, used to compute the 
phase diagram in Fig.[2]-a-b and[3] Here, hf, — sin(7Dc) is a single mode of the Laplacian operator. Thus the 
stream function y/ and the potential vorticity field q H are simply proportional to h/, for any values of the 
initial enstrophy Zo and the initial energy E. 





























































































38 



Fig. 3: Phase diagram of Fig.[2}a with hypothetical trajectories of the shallow water system in the presence 
of dissipation. Trajectory A-B: the dynamics dissipates small scale fluctuations of height and divergence 
field only. Trajectory A-C: the dynamics dissipates small scale fluctuations of potential vorticity. Trajectory 
A-D: the dynamics dissipates small scale fluctuations of all the fields. 


5.2 The effect of energy dissipation and enstrophy dissipation 

The actual shallow water dynamics is known to be characterized by shocks that prevent energy conserva¬ 
tion. In addition, the presence of viscosity, no matter how small it is, may lead to enstrophy dissipation, and 
more generally would break the conservation of potential vorticity moments. The aim of this subsection 
is to discuss qualitatively the effect of these dissipative processes on the large scale flow, assuming that 
the system evolves through a sequence a equilibrium states, which is a natural hypothesis if there exists a 
separation of time scales. 

For the sake of simplicity, let us focus on the phase diagram obtained in the energy-enstrophy ensem¬ 
ble and described in the previous subsection. Let us consider that the system has reached at some time an 
arbitrary equilibrium state denoted by A in Fig. [3] 

Let us first consider a case where energy is conserved, but enstrophy is dissipated. Then the system will 
evolve form point A to point C of Fig. 0 In other words, the small scale enstrophy Zt will be dissipated 
so that the enstrophy of the system tends to the minimum admissible value of enstrophy Z 2 m in(E m f). Note 
that the curve Z 2 min{E m f) is nothing but the curve of energy minima for a given Zt < Z/,. We saw previ¬ 
ously that the mean-flow energy E m f dominates the fluctuation energy Ef/ uc on this line. We conclude that 
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enstrophy dissipation alone drives the system towards a large scale flow without small scale fluctuations. 
This large scale flow vanished when topography vanished (since E nbx and Z b would also vanish). 

Let us now consider that the enstrophy is not dissipated and look the effect of energy dissipation. Let 
us first explain why one may expect that even weak dissipation can lead to a significant decrease of the 
large scale flow energy E m f when Efi uct ^ 0 (in the absence of small scale fluctuations, E m f would not 
decrease significantly with weak dissipation). At equilibrium, the energy of the small scale fluctuations 
Efiuct should be equipartitioned among all the modes of the height field and divergent field. Since there 
is an infinite number of such modes, this means a loss of energy through dissipative process, no matter 
how small they are. Since E m f decreases with Efi uct , dissipating the energy Eji uct amounts to diminish 
the energy E m f. For a given enstrophy Z 2 > Z b , as in the case of point A Fig. (Qjl, this dissipative process 
drives the system towards a zero energy state B. For a given enstrophy Z 2 < Z b , this process would drive 
the system towards the line of minimal energy. 

We expect to see both enstrophy and energy dissipation working together, so the trajectory of the 
system in phase diagram will be somewhere between the trajectory A — > B and the trajectory A —> C, e.g. 
the trajectory e.g. A —» D. We conclude that in the presence of topography and small scale dissipation, we 
eventually reach a geostrophic regime at large time, provided that the initial enstrophy is sufficiently low 
(otherwise the final state contains no large scale flow). In addition, each time the system is perturbed by 
adding a little amount of energy without changing the enstrophy, it drives the system a bit more towards 
the mixed state (E m f = E m j x . Z 2 — 0). 


5.3 Flow structure of the equilibrium states: application to the Zapiola Anticyclone. 

So far we have discussed energy partition for shallow water equilibria in the quasi-geostrophic limit. Flere 
we focus on the structure of the large scale flow associated with these equilibrium states. We consider for 
that purpose an oceanic application to the Zapiola anticyclone. 

The Zapiola anticyclone is a strong anticyclonic recirculation taking place in the Argentine basin above 
a sedimentary bump known as the Zapiola drift Ii54ll38ll24l . The anticyclone is characterized by a mass 
transport as large as any other major oceanic current such as the Gulf Stream. It is a quasi-barotropic 
(depth independent) flow, with typical velocities of the order of 0.1 m.s~ 1 , and a lateral extension of the 
order of 800 km. 

It is known from the earlier statistical mechanics studies that positive temperature states in the energy 
enstrophy ensemble of one layer models lead to anticyclonic circulations above topography anomalies, 
see e.g. [35] and references therein. A generalization of this result to the continuously stratified case with 
application to the Zapiola anticyclone is given in (481 . 

We have shown in this paper that quasi-geostrophic equilibria characterized by positive temperature 
states are also shallow water equilibria. The Zapiola anticyclone can therefore be interpreted as an equi¬ 
librium state of the shallow water model. Let us now show the qualitative difference between low energy 
states (when Z 2 > Z b and E m f —)■ 0) and high energy states (when E m t —> E mix ), which are shown on Fig. 
|4]-a and Fig|4]-b, respectively. 

In both cases the bottom topography is the same, and its isolines are visualized with thin black lines. 
Bottom topography has been obtained from data available online and described in B39II . We isolated the 
largest close contour defining the Zapiola drift (the sedimentary bump above which the recirculation takes 
place), and considered the actual bottom topography inside this contour, and a flat bottom outside this 
contour. It allows to focus on the interesting flow structure occurring above the Zapiola drift. 

Energy-enstrophy equilibrium states are then computed using Eq. (1234b in Appendix [F] We have seen 
that low energy limit corresponds to > +°° with Z 2 >Z b . Taking this limit in Eq. (1234b yields y/mf — 
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Fig. 4: Plots of the streamfunction isolines (colored lines) from higher values (red) to lower values (blue) 
over the Zapiola Drift topography iso-contours (black lines) for a small mean-flow energy (a)) and for a 
high mean-flow energy (b)). 


(7n — Zb) hb/[5 : the equilibrium state is a Fofonoff flow 1121 , meaning that streamlines are proportional to 
the isolines of topography, just as in Fig.[3}a. 

The maximum energy state corresponds to the case fi — 0, which corresponds to the mixed state 

defined in Eq. ( 1154b . The streamfunction is !//■„„> = (R 1 — A ) li/,. The operator (R 2 A) with R ~ 1 
is expected to smooth out the small scale topography features, just as in Fig. \4[a. 


6 Conclusion 

We have presented in this paper analytical computations of equilibrium states for the shallow water system, 
giving thus predictions for the energy partition into small scale fluctuations and large scale flow. Our results 
rely on the definition of a discrete version of the shallow water model. Once our semi-Lagrangian discrete 
model was introduced, the whole machinery of equilibrium statistical mechanics could be applied. 

We found that equilibrium states of the shallow water system are associated with the concomitant ex¬ 
istence of a large scale flow which is a stationary state of the shallow water dynamics, superimposed with 
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small scale fluctuations that may contain in some cases a substantial part of the total energy. The novelty 
of our work was to explicitly compute the contribution of these small scale fluctuations, and to decipher 
the physical consequences of the presence of these fluctuations. 

In particular, we found that the presence of small scale fluctuations implies a positive temperature for 
the equilibrium state. This explains a previous result by Warn II53II . who showed that equilibrium states in 
a weak flow limit admit no large scale flow when there is no bottom topography and no lateral boundaries. 
We have generalized these results, by showing that a large scale vortical flow exists at equilibrium when 
there is both rotation and bottom topography, or when there is a non-zero circulation. 

In the limit of weak height fluctuations, we found equipartition of the small scale kinetic and potential 
energy. We also obtained an interesting physical picture of the equilibrium state, which may be interpreted 
in that limit as two subsystems in thermal contact. One subsystem is the “large scale” potential vortical 
flow whose entropy is closely related to the one introduced heuristically by Chavanis and Sommeria (see 
|f7j). Our work provides therefore a microscopic justification of their entropy with a complete statistical 
mechanics derivation and a generalization of this results by including the presence of small scales fluctu¬ 
ations of the height and divergence fields. We note however that it is wrong to interpret the "large scale" 
potential vortical flow entropy as the entropy of the system, it is only one part of it. The other subsystem 
contains the field of height fluctuations associated with small scale potential energy, and the field of ve¬ 
locity fluctuations, associated with small scale kinetic energy. These velocity fluctuations are due solely to 
the divergent part of the velocity field. Warn obtained a similar result in the weak flow limit, by projecting 
the non-linear dynamics into eigenmodes of the linearized dynamics 1531 . He found an energy partition 
into a vortical flow on the one side, and on inertia-gravity waves on the other side, with a weak coupling 
between both subsystems. Hence we may interpret local height and divergent velocity fields fluctuations 
appearing in our model as inertia gravity waves. 

We studied the quasi-geostrophic limit for the large scale flow, taking into account the presence of 
small scale fluctuations. We recovered in this limit the variational problem of the Miller-Robert-Sommeria 
theory, with the additional constraint that the temperature is positive. We obtained phase diagrams in the 
particular case of energy-enstrophy equilibria, with explicit prediction for the ratio of energy between 
small scale fluctuations and a large scale quasi-geostrophic flow. This allowed to discuss the qualitative 
effect of small scale dissipation and shocks on the temporal evolution of the system. The main result is that 
such dissipative processes drive the system towards a minimum energy state (depending on the enstrophy), 
which may be non-zero. 

In view of those results, the semi-Lagrangian discrete model seems to be a good discretization of the 
shallow water system. It has the key desired properties to take into account of the conservation of fluid 
particle volume, while working in an Eulerian framework. This gives a clear framework, which allows for 
a rigorous derivation once the discretization is assumed. Moreover we stress that it leads to equilibrium 
states that are stationary states of the shallow water model, by contrast with other choices of discretization. 
What is not completely satisfactory however is that there is a degree of arbitrariness in the definition of 
the discrete model. We have tried to have the model consistent with the geometric constraints related to 
the Liouville theorem, even though the link between the model and the Liouville theorem is clearly not 
rigorous. There is clearly room for improvement, but we are afraid we are faced with extremely tricky 
mathematical problems. Nevertheless we guess that the invariant measure of the discrete model converge 
to the actual invariant measure of the shallow water system in the continuous limit, but proving this is 
beyond the scope of this paper. 

The statistical mechanics prediction of a vanishing large scale flow in the absence of boundary and 
bottom topography seems to contradict some numerical results performed in such configurations, in which 
case long lived vortex were reported, see e.g. ifTTIl . This issue is related to the estimation of a time scale 
for the convergence towards equilibrium. Indeed, if the coupling between the large scale flow (potential 
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vortical modes) and the small scale fluctuations (inertia-gravity waves) is weak, and if one starts from 
an initially balanced and unstable large scale flow, then this large scale flow may self-organize sponta¬ 
neously on a short time scale into an equilibrium state of the quasi-geostrophic subsystem. This justifies 
the physical interest of the variational principle for the large scale flow introduced by f7| without small 
scale fluctuation. According to our statistical mechanics predictions, the energy of this large scale flow 
should leak into small scale fluctuations, but this process may be slow if coupling between both subsys¬ 
tems is weak. This difficult issue was already raised by Warn |53j, and the interaction between geostrophic 
motion and inertia-gravity waves remains an active field of research BD- In particular, very interesting 
models of interactions between near-inertial waves and geostrophic motion have been proposed 15811141 
ED. In the context of statistical mechanics approaches, it has been proposed to compute equilibrium states 
with frozen degrees of freedom f36l or restricted partition functions in order to avoid the presence of 
inertia-gravity waves [Sl- 

One of the main interests of the present study is the prediction of a concomitant large scale energy 
transfer associated with a large scale potential vortical flow and a small scale transfer of energy that is lost 
into small scale fluctuations interpreted as inertia-gravity waves. In that respect, the shallow water model 
lies between three dimensional and two-dimensional turbulence. There are other models for which the 
energy may be partitioned into a large scale flow and small scale fluctuations. There are, for instance, some 
strong analogies with the case of three-dimensional axisymmetric Euler equations, that is also intermediate 
between 2D and 3D flows. A statistical mechanics theory has been recently derived for this system by 
Thalabard et al fidl . This statistical theory can also be understood as two subsystems in contact, one of 
them being the fluctuations, and thus leading to positive temperatures HD . From this key observation, HTj 
concluded that the temperature of the system, that measures the variance of the fluctuations, is positive 
and uniform in space, as previously stressed in 11271 . More generally, and beyond those simplified flow 
models that allow for analytical treatment, it is common to observe in bounded laboratory experiment the 
emergence of large scale structures at the domain scale superimposed with small scale energy fluctuations, 
see e.g. 1421 . It is not clear whether equilibrium theory may be relevant to describe such problems, but we 
believe that it is at least a useful first step to address these questions. 


A Invariant measure and formal Liouville theorem 

A. 1 Formal Fiouville theorem for the triplet of fields (h, hu. hv) 

The existence of a formal Liouville theorem for the shallow water dynamics is shown in this appendix. The shallow water 
system is fully described by the triplet of fields (h, hu, hv). We consider a measure written formally as 

d,u = CS> [h] & [hu] & [hv], 

with uniform density in (h, hu, hv)- space (C is a constant). The average of any functional A over this measure is 

Vs/[h,uh,vh], (s/) lx = j djU s /. (161) 

The term f 0 [h] & [hu] [hv] means that the integral is formally performed over each possible triplet of fields (h, hu, hv). 
The measure is said to be invariant if 

Vs/, 1(^/) M =0 (162) 

This yields the condition 

Vs/, J & [/i] 3> [hu] 9 [hv] J dx ~^d,h+ (hu) + d, (hu) = 0. 


(163) 
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An integration by parts yields 


'is/, J 3[h]3 [hu] 3 [/tv] srf J 


dx 


8d,h 89, (hu) 89, (/tv) 
Sh 8 (hu) 8 (hv) 


We say that the equation follows a formal Liouville theorem if we formally have 


/ 


dx 


f 89,h 89, (hu) 89, (hv) 

\ 8li + 8 (hu) 8 (hv) 


= 0 , 


which ensures that the measure dfi is invariant. 

The shallow water equations 0 and 0 can be written on the form 


= 0 . 


(164) 


(165) 


9, (Im) = —9 X ( hu 2 + -g/t 2 ) — dy ( huv ) + fhv, 


(166) 


We see that 


9, (hv) = —9 V (hv 2 + —gh 2 j — 9 X (huv) — fhu, 
9,li + 9 X (hu) + 9 y (hv) = 0. 


89,h 89, (hu) 89, (hv) 

~Sh~ - 1 " 4 ‘ 


= -V- 


5 (hu) 8 (him) 8(hv u) 

8h + 8 (hu) 8 (hv) 


(167) 

(168) 

(169) 


8 (hu) 8 (hv) 

As the divergence operator is a linear operator, it commutes with the functional derivatives. This allows to conclude that the 
measure /i is invariant. This shows formally the existence of a Liouville theorem for the fields (h, uh , vh). 


A.2 Change of variables from (hjiujiv) to (h. q.Li ) 

The microcanonical measure can formally be written 


dj Uh,hu,hv = @[h]@ [hu] 9 [hv] 8 {£ - E) 8 (J* - Z k ), 

k=0 

The constraints are more easily expressed in terms of the variables 


h, 


9 = 


-9yll+9 X V + f 
h 


^=A l ? 2 (9 x u + 9yv). 


(170) 


(171) 


It will thus be more convenient to use these fields as independent variables. We call J [(h,hu,hv) / (h,q, /x)] the Jacobian of 
the transformation. We proceed step by step to compute this Jacobian. The change of variables (h, hu, hv) —> (h, u, v) involves 
a upper-diagonal Jacobian matrix at each point r: 


J 


(h,hu,hv) 
(h, u, v) 


/ 1 u v \ 
0 h 0 
\00hJ 


which implies det(J[(h,hu,hv)/(h,u,v )]) = h 2 and 

& [h] & [hu] ^ [hv\ =h 2 S> \h\ 3) [«] 3) [v]. 


(172) 


(173) 


The change of variable (h,u,v) —» (/;, to,/r) involves linear operators that do not depend on space coordinates, thus the 
determinant of the Jacobian of the transformation is an unimportant constant: 


3> [h] 3) [u] 3 [v] = C3 [h] 3 [to] 3 [ju]. 


(174) 


Using a = qh— 1, the change of variable (/;, £0,/t) —J (h,q,[i) involves an upper diagonal Jacobian matrix at each point r: 
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J 


'{h,co,nY 

{h,q,n) 


1 q 0 
0 h 0 
0 0 1 


(175) 


with a determinant det (J[(h, / (q,h,n)\) = h. Finally, the Jacobian of the transformation is J [(hu.hv.h) / (q,h,n)\ = /t 3 

and the microcanonical measure can formally be written 


= Ch 3 ® [A] ® [<?] # [At] 5 - E ) fl 5 (3* - Z,). (176) 

k=0 

We note the presence of the pre-factor /; 3 which gives the weight of each microscopic configuration in the space. 


B Relevance of the constraints for the discrete model 

In this appendix, we explain how the dynamical invariants of the shallow water model, given in Eqs. 1321 and 1341 respec¬ 
tively, are related to the constraints of the microcanonical ensemble for the discrete model, given in Eqs. l92t and l97l i. 
respectively. 


B. 1 Areal coarse-graining for continuous fields 

Let us con sider a field g (x) on the domain & where the flow takes place, and let us consider the uniform grid introduced in 
subsection l3.2.1l We define the local areal coarse-graining of the continuous field g (x) over a site (;'. j) as 

8ii = N f dxg(x), (177) 

J site (;j) 

where 1 /N is the area of the site (i, j) and where J s j te ., means that we integrate over the site (i, j) only. With an abuse of 
notation, we use here the same notation as in Eq. 1501 . since the coarse-graining operator defined in Eq. 11771 generalizes 
to the continuous case the areal coarse-graining operator defined in Eq. {50) for the discrete microscopic model, taking into 
account the fact that for a fluid particle “n” of area dx„ and height h n , we get N dx„ = 1/ (Mh„). 

We denote g the continuous limit (large N) of g„ . Integrating a continuous field g amounts to perform the integration 
over its local average field g: 

J dxg(x) = f dxg(x). (178) 


B.2 Potential vorticity moments 

Using 11781 . the potential vorticity moments in Eq. 1341 simply leads to 


2 * = 



(179) 


Now that the potential vorticity moments are expressed in terms of the areal coarse-graining of moments of h and q, it can 
directly be expressed in terms of the probability density field p (<j/,, a q , ay,) through Eq. 1901 . and we recover the expression 
of the constraint given in Eq. |92), whose discrete representation is given in Eq. 1561 . 
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B.3 Energy 

Using GUl i, recalling that we restrict ourself to bottom topographies such that hg = hg, the total energy of the shallow water 
model defined in Eq. {32) can be decomposed into a mean flow kinetic energy defined in Eq. ( I94t . a potential energy term 
due to local height fluctuations and defined in Eq. 1961 . a fluctuating kinetic energy term 

£’cjiuci = ^ ~ S m f - S 5h (180) 

Scjluct = \ J dx (/iu 2 -Tm; nf ) , (181) 

where the velocity fields u and u„,f are computed from the triplet (h,q,fx) and from the triplet ( h,hq,h)i ), respectively 
through 

u = V-V + V<)>, [hq — f) = Ay, n = A l / 2 <j), (182) 

and through 

u »,/ = V (hcj /( = A \ffmf, ~=~ A / (j) m f . (183) 

h 

We want to discuss the relation between decomposition of the energy for the discrete model, Eq. {97j and the decomposition 
for the actual total energy defined in Eq. <321 . Our construction is relevant if these two decomposition coincide in the 
continuous limit, or equivalently if S c ji uct is equal to {95} in the continuous limit. In the following we show that this is 
the case if some cross correlations are actually negligible. More precisely, we assume that 

1. For any p ositive integers k,l,m , the coarse-grained fields h k q l defined through the c oarse graining procedure in 
Eq. ( 11781 ) exist. In the framework of our microscopic model introduced in subsection l3.2l this hypothesis is automat¬ 
ically satisfied by assuming that the cut-off fl m in,Hmax,<lmm,qmax,hmax scales as N a with a < 1. Other fields may be 
characterized by local extreme values such that the limit defined in Eq. ( 1178b does not converge. For instance, we will 
see that the actual divergence £ = A</> of the equilibrium state is not bounded, i.e. that £ would have no meaning. 

2. The fields h, q, / 1 are decorrelated (in particular, h k q' n' n = /t^'/T"). This point will be shown to be self-consistent when 
computing the equilibrium state. 

3. The coarse-grained divergent velocity field is equal to the mean-flow velocity field V<j> = 

4. A | [tj) — 0 m y) 2 j = 0. While V(0 — (j) m f) is a random vector field characterized by wild local fluctuation, this hypothesis 
amounts to assume that those fluctuations have no preferential direction. 

We believe that these four assumptions would be satisfied by a typical triplet of fields ( h,q,fx ) picked at random among all 
the possible states satisfying the constraints of the dynamics. By typical, we mean that an overwhelming number of fields 
would share these properties. 

We then prove that these four assumptions are sufficient to prove that <S c ji uct is equal to S’g^ ( 195b in the continuous 
limit. According to the assumption 1, W = hq— 1 is well defined. Classical arguments show that the streamfunction of the 
coarse-grained vorticity field to is equal to the streamfunction of the vorticity field, i.e. thatVl jf m f = Vi/t, see e.g. mm. 
Qualitatively, this is due to the fact that inverting the Laplacian operator smooth out local fluctuations of the vorticitjQ This 
yields 


U = U m f + V (0 — . ( 184 ) 

Injecting this expression in the kinetic energy density expression Eq. {HD- using that h and fi are not correlated (assumption 
2), and using V0 = V<j> m f (assumption 3) yields 


Scjluct = dx/?(V (0 - (j> mf )) 2 . 


(185) 


13 The divergent part of the velocity field can not be treated in the same way. Indeed, the operator A x ! 2 is less smooth 
than the operator A ~ 1 , and one can not derive <j> = <j) m f by inverting /i = A 1 Z 2 ij> . One may want to use f = A (j> , but the result 
would be the same since £ is not bounded (the field fi is characterized by fluctuations which are controlled by the kinetic 
energy, and hence by the total energy, but this is not the case for Q. 










46 


Let us now remember the definition of the coarse-graining operator in Eq. OH): 

(V (0 - 0mf)) 2 = J dx (V (0 - 0mf)) 2 , for x e Sij, (186) 

Sfj 

where Sjj is the surface covered by a grid site j). An integration by parts yields 

N ( dx (V ( 0-0m f )) 2 = —N f dx [(a~ 1/2 (h -ju)) (A l/2 (p -p))] +N (f> dl n- (0-0 mf ) V (0-0 mf ) (187) 

Sij s u 3S,j 

Proj ecting the first term of the rhs on Laplacian eigenmodes allows to simplify the expression of the first term of the rhs in 
Eq. GO': 


— N J dx [(zl 1/2 (p-p)) (zl 1/2 (p-p))] = Njdx (p-p) 2 . (188) 

Sij Sij 

The second term of the rhs in Eq. d 1 87t can be written as 

N <£ dl n- (0 - 0 mf ) V(0- 0mf) 

3 Sij 

which, according to assumption 4, vanishes in the large N limit, 
simply expressed as 

Scjiuct = \ Jdxh (V -]u 2 ) . (190) 

Finally, we use again the assumption 2 to get 


= ^ f A[(0-0 mf ) 2 ] (189) 

J Sjj 

Finally, the kinetic energy density of the fluctuations is 




C,fillet ' 





(191) 


which is the expected result. 


C Critical points of the mean-flow variational problem 

In this Appendix, we compute the critical points of the mean-flow variational problem i ll 001 stated in section[3] In a first step, 
we solve an intermediate variational problem in order to show the factorization of the probability density p with a Gaussian 
behavior for the divergence fluctuations. Knowing that, we solve in a second step the original variational problem. 


C.l Intermediate variational problem 


As the energy and the potential vorticity moments depend only on the coarse-grained fields h. h 2 , hji. h[i 2 and the local 
potential vorticity moments j, we introduce an intermediate variational problem where these coarse-grained fields are 
given as constraint: 


max 

pJp = 1 


y[p\ 


h,h 2 ,hn,hn 2 ,^hq k ^ "fixed" 


(192) 


The idea of introducing the intermediate variational problem is to find a sim pler ansatz for the probability density field p . 
This ansatz will be used afterward into the general variational problem ( 1 1 001 ). 
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In order to compute the critical points of the variational problem {J92J, we introduce the Lagrange multipliers a* (x), Cfe (x), a/ Ifi (x), ( 4 ^ 2 (x), {( 4 ^ (s 
associated with the constraints h , h 1 , hji, hp 1 , ^hq k | and the normalization constraint, respectively. Using Eq. jS9| and 


the first variations 

V5p, 55- / dx 

leads to 


/‘ 


a h 8h + a h2 8h 2 + a hp 8hn + a hfl2 8hfi 2 + V a hq ^8hq k 

k= 1 




8p 


= 0, 


(193) 


CT/,log ( — ) + o h +t, + a h a h + a h2 al + a hll a h cfy + ai^OhOjl + £ a hq ,kOhO q = 0 


(194) 


We readily see from Eq. i ll 94) that the probability density p factorizes into three decoupled probability densities p q , pi, and p^ 
corresponding respectively to the probability densities of the potential vorticity, the height and the divergence: 


p = Pq(x,o q )p h {x, a h ) p M (x, cfy). 


(195) 


Using the constraints p = f defy (fyp^ , p 2 = / defy efyp^, , as well as the normalization constraints f defy p^ = f defy pi, = 
/ d <fy Pq = 1, we get 


Pq ( x i <fy) — 


exp -£ a hq , k (x)o k q 


J d<fy exp a hq .k (x) <fy‘ 


P/,(x,ct/,) = 


er ; ‘exp -a h2 (x)a h - 




/ d °/; < 


exp -a /,2 (x) a; - 


<Jh 


(196) 


exp - 


P/i ( x , op) — 


1 (fy-p(x)) 2 


2 p2(x)-p 2 (x) 


(2fy 1/2 (p 2 (x)-p 2 (x)) 


T72 


We could now re-inject these expressions into the main variational problem ( 1 1 00b . but only factorization and the Gaussian 
form of p t , will be kept as an ansatz for p, which will simplify the computations. Thanks to this intermediate variational 
problem, we now know that the critical points of the original variational problem must be of the form: 


P ( x ,C4,(fy,(fy) = Ph (o), , x) p q (<fy, x) 


exp 


1 (ty-Ap 2 


•_ x 1/2 • 

.. 0 —2 \ 


(2tt) 1/2 (p 2 -P 2 ) 


(197) 


The entropy defined in Eq. ( f89l> is therefore (up to a constant): 


y [p/np?.p.p 2 -p 2 ] =- f dxdcj/, a,,pi, log — j dx/? J dtJqPq log (p^) T J dx ^ log (p 2 - p 2 ) . (198) 


As a consequence of Eq. fT97t . the height field, the potential vorticity field and the divergence field are decorrelated. 
This property allows to rewrite the energy defined in Eq. d93t 


[p/,,p, 7 ,p,p 2 -p 2 ] =g mf \h,q,p] +S Sfl [/qp 2 -p 2 ] +£ Sh [/i,/i 2 ] 


(199) 
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where 

<4/| Ti.Ji.~p] = 3 /dx (Tm; nf +g(h + h b -iy^ 

< S 5n JX /x 2 - p 2 ] = \ f dxh (p 2 - p 2 ') (200) 

#sh yh h 2 ]= if dx (h 2 - Tfj 

with U;„/ = V-'-zW 1 fqh — /) + V4 ~ 1/<2 ]u . Similarly the potential vorticity moments ((96} can be rewritten 

Vk J21 [S,?] = J dx hf (201) 

where the coarse-grained moments are now defined as 

h l = J do;, aj,p h , q 1 * = J Aa q a"'p q . (202) 

Thus, the general variational problem of the equilibrium theory given in Eq. Hoot can be recast into a new variational 
problem on the independent variables pi, (x, O/,) , p q (x, a q ) , p (x) and |X 2 — p 2 J (x): 

S(£,D)= max _ jj?' Ip ; ,,p, ; ,p,p 2 - p 2 ] | S [p,„p 9 ,p,p 2 - p 2 ] = £, V* J* [p h ,p q ] = Z k J • (203) 

Pl,,Pc,,HS 2 -H 2 
I Ph='j P?=> 


C.2 Computation of the critical points 


In this subsection, we compute the critical points of the variational problem defined in Eq. 12031 . We introduce the Lagrange 
multiplier /3, {a,t}j. >0 , %q{ r ) and <*/, (r) associated respectively wi th th e energy, the potential vorticity moments and the 
normalization constraints. Critical points of the variational problem 12031 are solutions of 

\J8p q ,8p h ,8p,8 (p 2 -p 2 ) . Sy-pSS’- £ a k 82? k - j Ax (^ q J 8p q + fy, J 8p = 0 . (204) 

The first variations of the macrostate entropy y ( fl98t are 


” sy 

^Ph 




sy 

8y 

W 


sy 

sy-P 2 ) 


= -a h (tog ^ J + 1 j - a h f d a q p q log {p q ) + a h \ log (p 2 - p 2 ^j 

= -Mi°g(p«) + i) 

= 0 

_ h 

2 (?-M 2 ) 


(205) 


First variations of the energy given in Eqs. 11991 and 12001 contain three contributions: 8S = 8S m f + 8Sg„ + 8£$ h . The 
first contribution is 

(206) 


8f&mf — j dx [B m f8h -j- (hu mf J • <5u„,/], 


where B m f = u 2 , ^/2 + g (/? + hg — l) is the mean-flow Bernoulli function defined in Eq. 11151 . Then, using the Helmholtz 

decompositions u m f = V 2 - iff,,, / + V </) m / and recalling that hu m / = V 1 - + V <P m f , two integrations by parts with the imper¬ 

meability boundary condition yield 


8&mf — j dx \B m f8h *Pmf 8A Ijt/ri f ( P,„ f8A0,„ f . 

Using A Ymf = hq — f and A l ^ 2 <j> m f = p and the definition of the operator A */ 2 leads to the final expression 

8S m f = J dx ^(B mf -q'Bjnf) 8h-h'F mf 8q-A 1 / 2 <l> mf 8p^ . 


(207) 


( 208 ) 
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Finally, we get: 


bs 

Spi, 

s,g 


3% 

bS 


W 


bS 




= o„ (fi m/ - q¥ mf + ^ + g (f - h )) 
Oq^mfh 
= -A l / 2 <P mf 

_ h 
2 


First variations of the potential vorticity moments are 


VAeN 


bsr k 

Sph 

b3fk 

8p« 

s -?k 

W 

SSi 

s(p 2 -p 2 


= O h q k 
= ho* 
= 0 
= 0 


(209) 


( 210 ) 


Injecting Eqs. <2051 . ( I209K and <21 01 1 in Eq. <2041 . and collecting the term in factor of 5 u i 2 — fi 2 ) leads to 


n 2 -~P 2 = 


P 


Injecting Eq. <21 11 1 in the expression of p^ given in Eq. <1961 yields then 


= ex p(-^ (a M -p(x)) 2 ^. 

Similarly, collecting the term in factor of Sp q in Eq. <2041 leads to 


0 = h (log(p ? ) + 1 ) -Po q %,fh+ 52 a k ho* + £, q 


which, using the normalization constraint, leads to 


p q (x, o q ) = — exp ( pOq¥ mf - £ a k o; 


G q = j Ao q exp o q % lf - 52 a k o^j . 


( 211 ) 


( 212 ) 


(213) 


(214) 


Note that the sum inside the exponential is performed from A = 1 to A = +°°. The Lagrange parameter t, q has been determined 
using the normalization condition for the pdf. 

Collecting the term in factor of 8p/, in Eq. <2041 yields 


- ( l0g (ff^) + *) “ / dcJ 9 p ? log ( p «) - ^ lo S(P)~P (^mf+^2~+g{y—h^ - + Pq% t f ~ Yk a kq k = 0, 

(215) 

which, using Eq. <2l4t . leads to 

~ ^l°g ^ +l°g(G ? ) _ 2 log ^ ^W + - Jr 1 ~ «o = 0. (216) 

Usi ng the fact that G ? and B m f are fields depending only on x, and using the normalization constraint for the pdf p/,(x, a),), 
Eq. ( 12 1 61 yields 

p h (x,o h ) = ^exp(-pjO h -^, G k = J d<7 ; ; Oft exp ■ ( 217 ) 

Injecting Eq. <2 1 7b back into Eq. <2 1 6b yields 

B m f = fi 1 log(G ? G*) +gh + p 1 Oq — 4. + -log/3^ . (218) 
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One can notice that cto the Lagrange parameter related to the conservation of the total mass appears only here. Thus the last 
term /3 1 (cto — 3/2 + log (/3) /2) in Eq. <2181 can be computed from the conservation of the total mass 5o = Zo and will be 
denoted Ao in the following. 

Collecting the terms in factor of 5JT in Eq. <204 l i leads to 

<*>»,/= 0. (219) 


D Global maximizers of the entropy of the large scale flow 

We compute in this appendix an upper-bound for the macrostate entropy of the large scale flow defined in Eq. fl30t . for 
a given set of potential vorticity moment constraints defined in Eq. <921 (and arbitrary energy), and then show that when 
Zi = / and lib = 0, this upper bound for the macroscopic entropy is reached by the rest state. 

This upper bound is the solution of the following variational problem: 

Smf ,max — max {y mf \h,p q ] | Vfc J* \h,p q ] =Z k }. (220) 

H,Pq 
fPq= 1 

Introducing Lagrange parameters {)^} < .>o associated with the potential vorticity moment constraints and the Lagrange pa¬ 
rameter t, (x) associated with the normalization constraint, the cancellation of first variations yields 

V8p q ,8h, 8y mf -Y,r k 83f k + fdx|5T=0. (221) 

k=0 J 


The solution of this equation is 


P‘, 


exp E*=i S 5 


— Pglobal ( & ) 


f drr exp E*=i 

where p g , 0 bai depends only on the potential vorticity moments constraints {Z k } lc>1 , and is independent from x and 


( 222 ) 


h (x) = 1 


(223) 


Note that the states characterized p q = p g i 0 b a u h = 1 are solutions of the variational proble m in E q. <2201 . but this is 
only a subclass of the solutions of the variational problem of the equilibrium theory given in Eq. (S3) , which includes an 
additional energy constraint 

We have shown in subsection (4 A) that for a given p 9 , the large scale flow which is a solution <1531 is obtained by solving 
Eqs. <1171 and <1181 for %,f and h. Here we consider the particular case p q = p g , 0 bai ar| d h= 1 . One can compute hq global = 
f dx apgiob ai = Z|. We conclude that the large scale flow of the equilibrium state is also a global entropy maximizer, i.e. a 
solution of <2201 when 


Z l -f = AW mf , (224) 

J(Vt/) 2 + s/i t =A 2 . (225) 

where = /3 logG 9 — A \ is a constant. We see that in the case (Z\ = /, h b = 0) , the solution of Eqs. <2241 and <2251 is 
the rest state %,f = cst (with An = 0). We conclude that the maximum of the macroscopic entropy of the large scale flow is 
reached by a flow at rest when there is no circulation (Zj = /) and no bottom topography ( h b = 0). 
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E Comparison with a Eulerian discrete model 

The aim of this appendix is to discuss the construction of a possible invariant measure for the shallow water equations 
through an Eulerian discretization. We prove that the obtained equilibrium states differ from the one obtained through the 
semi-Lagrangian discretization used in the core of the paper. Moreover, we prove that the equilibrium states are not stationary 
states of the shallow water equations and that the statistical equilibria are not stable through coarse-graining. 

We define a purely Eulerian discrete model by considering the same uniform N x N grid as for the semi-Lagrangian 
model, but assuming that each node is now divided into a finer n x n uniform microscopic grid. A microscopic configuration 
is given by the values of the fields ( h , q, fi) for all the nodes of the microscopic grid: 


Vmicro — J .ij • d t J.ij •. PU.ij } 1 (226) 

1 <i,j<n 

where ( I.J ) and (i.j) correspond respectively to the position on the macroscopic grid and the position on the microscopic 
grid within the macroscopic node. 

Contrary to the semi-Lagrangian model, the Eulerian model has the desired property to possibly be compatible with 
the formal Liouville theorem derived in Appendix [A] for the continuous dynamics (although no mathematical result exist). 
However, the volume of fluid varies from one microscopic grid node to another in the Eulerian model, depending on the 
value of the height By comparison, our semi-Lagrangian approach respects the Lagrangian conservation laws (the 
height h is defined through the particle mass conservation). Because of the need to go through a discretization to build the 
microcanonical measure, we see that both the Eulerian and the semi-Lagrangian approaches necessarily break part of the 
geometric conservation laws of the continuous model. Hopefully rigorous mathematical proof of the convergence of the 
measures of one of the discretized model to an invariant measure of the continuous equations will settle rigorously this issue 
in a near future, however nobody seem to know how to attack this problem mathematically. We are thus led to the conclusion 
that based on current knowledge, there is no clear mathematical or theoretical a priori argument to choose either the Eulerian 
or the semi-Lagrangian discretization in order to guess the microcanonical measures. For now, the use of one discrete model 
or another to guess the microcanonical measure of the continuous shallow water equations can therefore only be justified a 
posteriori. 

Let us now define the empirical density field as 

1 " 

) ty m icro] — (Jt/jjj rT/,) <5 {q IJij (Jq ) S (p/jjj Oji). (227) 

n i,j= 1 

One can now compute the entropy of the macrostates p = {y m i cro | W, J du [y m , cm ] = Pu}, which, after taking first the limit 
n —>■ oo and then the limit N —» °° leads to 

^Eul [p] = - J dxd(7;,dcr f; dO)j p (x, 0 ),, O q , <T M ) log ^ P( X ’ °>»Oq, ^ ( 228 ) 

This Eulerian macrostate entropy has to be compared with the macrostate entropy for the semi-Lagrangian discrete 
model given in Eq. <89b . We can switch from expression to the other by changing p into CT;,p. We note that the two entropies 
become equivalent at lowest order in the limit of weak height fluctuations and weak height variations. However, in the 
general case, they are different, and therefore lead to different equilibrium states. In particular, is it straig htforward to show 
that because of the absence of the factor C/, in the expression of this Eulerian macrostate entropy j228b . the critical points 
p(x,t7/,,(7 9 ,CTp) of the microcanonical variational problem do not factorize. Consequently, small scale height and velocity 
fluctuations of the equilibrium state are correlated. One can then show that those correlations are associated with non-zero 
Reynolds stresses in the momentum equations. In particular, the equilibrium state of Eulerian model satisfies 


J('P,q)=-J(W',q')-J(4>\q'), (229) 

where the r.h.s. is non-zero. If one removes those small scale fluctuations, the large scale flow is not a stationary state of the 
dynamics since 7( l f',g) ^ 0. In other words, the equilibrium states of th e Eul erian model are not stable by coarse-graining, 
contrary to t he eq uilibria of the semi-Lagrangian model. Moreover, Eq. ( 12291 ) and the properties of stationary states derived 
in subsection l2.2[ implv that neither the potenti al vorticy field q nor the Bernoulli potential B m f can simply be expressed as a 
function of As shown in subsection l4.3.1b . the fact that B m f is a function of is essential to prove that the equilibrium 
is characterized by geostrophic balance at lowest order in the Rossby number Ro, when Ro —> 0. Consequently, the proof of 
geostrophic balance derived in the framework of the semi-Lagrangian model does not hold in the framework of the Eulerian 
model, unless the bottom topography is sufficiently small (h/, ~ Ro). 
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Let us finally argue that the stability by coarse-graining is a desirable physical property for the equilibrium states. 

The first argument is a body of empirical observations. In either experiments, geophysical flows or numerical simulations 
flows governed by the shallow water equations (or the Navier-Stokes equations or the primitive equations in a shallow water 
regime) in the inertial limit (when they are subjected to weak forcing and dissipation, with a clear time scale separation) 
do actually self-organize and form large scale coherent structures for which there is a gradual decoupling of the flow large 
scales and small scales. A prominent example is the velocity field of Jupiter’s troposphere. 

The second argument follows. Macrostates that evolve through an autonomous equation, must increase the Boltzmann 
entropy. This is a general result in statistical mechanics, which is a consequence of the definition of the macrostate entropy 
as a Boltzmann entropy. Indeed as the Boltzmann's entropy measure the number of microstates corresponding to a given 
macrostate, it must increase for most of initial conditions. When there is furthermore a concentration property (which is the 
case for the shallow water case, both the the Eulerian and sem-Lagrangian discretizations), the number of initial conditions 
for which the entropy can decrease decays exponentially with N (N is often the number of particles in statistical mechanics, 
here the number of degrees of freedom of our discretization). As a consequence, the set of equilibrium macrostates (entropy 
maxima) has to be stable through the dynamics for most initial conditions. In the shallow water case, in statistical equilib¬ 
rium, obtained either using the semi-Lagrangian or Eulerian discretization discussed above, the stream function concentrates 
close to a single field (the stream function fluctuations vanish in the large N limit). As a consequence the macrostate stream 
function, which is a single field thanks to this concentration property, has to be stationary for the dynamics. Those two prop¬ 
erties, that follow from the definition of the Boltzmann entropy, are actually verified for the equilibrium measure constructed 
from a semi-Lagrangian discretization, but not for the equilibrium measure constructed from a Eulerian discretization. For 
this reason, we conclude that the microcanonical measure constructed from the purely Eulerian discretization is inconsistent 
with the shallow water dynamics. 

We note moreover that the stability of the equilibrium macrostates through coarse graining ensures that the equilibrium 
states of the inviscid system are not affected by perturbations such as a weak small scale dissipation in momentum equations. 
This property is not a-priori required for the invariant measure of the shallow-water equations. However It is extremely 
interesting as it is a hint that this invariant measure may be relevant for non perfect flow in the inertial limit. 


F Energy-Enstrophy ensemble 

F. 1 Computation of the critical points 

In this Appendix, we compute the solutio ns of the variational problem CJD and describe the corresponding phase diagram. 
Critical points of the variational problem ( 1158b are computed through the variational principle: 

V5p ? , 8y’ mf , g --^—8g mf , g -y 2 82f g2 -y l 82f gl - [ dxt;(x) [do q 8p g = 0, (230) 

Zfluct J J 

where yi, y\ and t, (x are Lagrange multipliers associated with the enstrophy conservation, the circulation conservation, and 
the normalization respectively. Anticipating the coupling between the large scale quasi-geostrophic flow and the small scale 
fluctuations, the temperature is denoted Efi uc , (the inverse temperature is the Lagrange parameter associated with energy 
conservation).This yields 


Pg ( x ’ °?) ~ 


2 7i (Z 2 — Z 2 ) 


ex P 


2(Z 2 -Z 2 ) 


=t- o a - 


y^mf 

Efluct 


-tt)(z 2 —z 2 )) 


where we have introduced the enstrophy of the coarse-grained potential vorticity 


Z 2 = 


/ 


dxql 


(231) 


(232) 


Injecting ( 123 1 1 in Eq. J 147b , using the mass conservation constraint given in Eq. i ll44b and the zero circulation constraint 
% [g ? ] = 0 yields 


q g = PYmf with/3 


(Z 2 -Z 2 ) 

Efluct 


(233) 


Note that /3 is necessarily positive given that Z 2 — Z 2 > 0. Injecting Eq. ( 1233b in Eq. ( 1149b . the streamfunction can be 
computed explicitly by solving 











53 


P tymf — ^ tymf Yinf hb- (234) 

In order to solve this equation, it is convenient to introduce the Laplacian eigenmodes of the domain with k G N + : 

Ae/i = —Xleii with e/< = 0 on dS>, (235) 


where the eigenvalues —A^ are arranged in decreasing order. We assume those eigenvalues are pairwise distinct. We also 
assume that the bottom topography is sufficiently smooth such that £ t . |/zm| 2 A^ < +°°. Then, given that /3 > 0, the projection 
of the mean flow streamfunction on the Laplacian eigenmode ejr(x) is obtained directly from Eq. d234t : 


Yk = 


hbk 


p+tf + R- 2 


(236) 


We see that there is a unique solution l j/ m f for each value of p. This solution is therefore the equilibrium state. All the 
large scale flows associated with statistical equilibrium states of the shallow water system restricted to the energy-enstrophy 

ensemble with zero circulation are obtained from Eq. 12361 when varying P from 0 to +°°. 


F.2 Construction of the phase diagram 


The problem is now to find which equilibrium state is associated with the constraints ( E,Z 2 ). In the following, we explain 
how to find the equilibrium states associated with parameters ( E m f,Z 2 ), and how to compute the temperature Ef/ uc for each 
of those states. It is then straightforward to obtained the total energy E = E m f + Efi uc . 

Injecting Eq. J236I) in the expression of the quasi-geostrophic mean-flow energy defined in Eq. fi50) yields 


En,f - 7 £ (^t 


+ R- 


I \hbk\ ^ 


(237) 


The mixing energy defined in Eq. d 155b is recovered for jS = 0, given that Z\ = 0. In the range /3 > 0, the energy E m f 
is a decreasing function of /3, varying from E m f = E nl i X to E m f — 0, see Fig. [5]-b and Fig.[5}d. 

Injecting Eq. d236b in the expression of the macroscopic enstrophy given in Eq d232b yields 


= £ \h bk \ 2 

/t=i 


1 -- 


(238) 


The potential enstrophy Zj, defined in Eq. is recovered for /3 = +°°. The macroscopic enstrophy Z 2 is an increasing 
function of P, varying from Zi=0 (for /3 = 0) to Z-± = Z* for (P = +°° ), see Fig. [5]-a and Fig.^Jc. 


Two expressions of the macroscopic enstrophy Z 2 in terms of the parameters P have been obtained: one is given by Eq. 
d238b . the other arises from the definition of P in Eq. {233), which yields 


Z 2 = Z 2 - Efi uct p. 


(239) 


For given values of Efi uct and Z 2 , the values of P and Z 2 are obtained by finding the intersection between the two curves 

defined in Eq. {238) and {239), respectively. Once p is obtained, Eq. d237b gives directly the value of the mean-flow energy 
E m f. The phase diagram presented in Fig.[2]is obtained numerically by using this procedure. Graphical arguments presented 
in the following allow to understand the structure of this phase diagram. 
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Fig. 5: a) Variation of the macroscopic enstrophy Z 2 over /3, case Z 2 > Z/,. b) Variation of the mean-flow 
energy E m r with )3, case Z 2 < Z fo . c) Variation of the macroscopic enstrophy Z 2 with /3, case Z 2 < Zfc. d) 
Variation of the mean-flow energy E m f with /3, case Z 2 < Z/,. 


F.3 Limit cases for the energy partition 

Let us first note through Figs. [3Ja and[5}c that /3 is an decreasing function of Efi uc/ . Indeed, jS is given by the intersection 
between the solid curve representing the expression of Z 2 given Eq. ( 1238b and the dashed line representing the affine ex¬ 
pression of Z 2 given Eq. d239b where —Efi uct is the slope. Then we know that the total energy E = E m f([ 5 ) +Efi uct is an 
increasing function of Efj uct . Let us now consider different limit cases. 

The limit E —» °° with Z 2 fixed: In this limit, we have also Ef !ucl —» and /3 — > 0 Hence, one gets from Eq. d237b (see 

also Figs. [5}b and[5}d): 

lim E mf =E mix , lim %^=0. (240) 

£->+°o ' £^+00 E 


The lowest E limit with Z 2 < Z|, fixed: In this limit, we have also Efi uct —> 0. One gets from Fig.|5]-c that /3 —» j5 nmx (Z 2 ). 


Hence, E, 


'. m f reaches a minimum admissible energy E m j n (Z 2 ) = E m f \J5 ma x (Z 2 ) j . Then: 


lim 


Emf _ 1 


E—>E m in(7 / l) E 


(241) 
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The limit E —» 0 with Z 2 > Z|, fixed: In this limit, we have also Ef\ uc , —)• 0. One gets from Fig.[5}a that j5 —> 00 and that 
Z 2 —> Z b . Hence, from Eqs. 1237b and 12391 . we obtain: 


Vmf = C b p- 2 +o(p- 2 ) 1 

v / with r,. = - 


l* £ ~, <Z! - Z ‘> £ 


/, with C 6 = -£|/t M |"(\ 2 +^ 2 ). 


(242) 


fluct 


Here, CJ, is a constant depending on the topography only. Thus we have E m f ~ El luct C b / {Z 2 — Z b ) 2 , which leads to: 


lim ■ 
£-> 0 £ 


"m/ 

fluct 


= 0 , lim —£ = 0 . 
£—>0 E 


(243) 


The limit E —> 0 with Z 2 — Zb ~ C a E“ with a > 0: In this limit, we have Efi uc , —> 0. One gets from Fig. [5}a that 
/3 —f 00 . Hence, from Eqs. 1237b and 1238b . we obtain: 


mf = C„)3- 2 + 0 (/3- 2 ) 
: b -Z 2 = 4C b j3- 1 +o(ji- 1 ') 


\E, 

Z h - 


(244) 


where C b is defined in Eq. 1242b . From those two equations along with Eq. 1239b and using Z 2 — Z b ~ C a E a , we can extract: 

C a E a + \J C^E 2a + 20C b E 


P 


2 E 


(245) 


Now, we have to consider different cases for the value of a. 

For a > 1/2, we have from Eq. 1245b that /3 ~ \Z5C b E~ 1 / 2 . Injecting this in Eq. 1244b . we gets: 

.. E m f 1 
lim —- = - 
£-> 0 E 5 

For a < 1/2, we have from Eq. 1245b that /3 ~ C a E“ _1 . Injecting this in Eq. 1244b . we gets: 

.. E m f 
lim —- = 0 

e -¥0 E 


(246) 


(247) 


For a = 1/2, we have from Eq. 1245b that /3 ~ C t j 2 E '/ 2 ^1 + y/l + 20C b /C]j l J /2. Injecting this in Eq. 1244b . we gets: 

2 C b /C 2 a 


lim — 

E-y 0 E 


(i+vi+20 c b /c*y 

Contrary to the previous cases, here, the partition of the energy depends on the bottom topography. 


(248) 
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